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This  study  was  an  evaluation  of  a  two-cloud  layer  heating 
model  for  inclusion  in  a  numerical  weather  prediction  model 
over  the  ice-free  ocean.   The  model  employs  empirical  formu- 
lations for  atmospheric  absorptivities ,  scattering,  and  re- 
flectivity, and  cloud  and  earth  surface  reflectivity  to 
compute  solar  insolation  received  and  absorbed  at  earth  and 
in  the  key  tropospheric  layers  above.   Long-wave  emissivity 
formulas  by  Sasamori  for  water-vapor  and  carbon  dioxide  are 
computed  to  determine  cooling  by  the  earth  and  troposphere. 
To  complete  the  heat  budget  scheme  for  the  earth's  surface 
and  near  surface  boundary  layer,  vertical  transport  of  sen- 
sible and  latent  heat  at  the  surface  was  evaluated. 

To  check  the  validity  of  this  model,  both  instantaneous 
heating  rates  and  24-hour  averaged  rates  were  computed  for 
the  earth-atmosphere  system,  the  atmospheric  column  and  the 
earth's  surface  using  data  in  the  form  of  grid  point  sound- 
ings over  the  Atlantic  and  Pacific  Oceans  for  15-16  October 
1973,  provided  by  Fleet  Numerical  Weather  Central. 


EY. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 16 

II.  DATA  PREPARATION 20 

A.  TEMPERATURE  TERMS ■, 24 

B.  MOISTURE  TREATMENT  IN  THE  SOUNDING 24 

C.  PRESSURE-SCALED  ABSORBER  MASSES  26 

D.  CLOUD  PARAMETERIZATION  28 

III.  TERRESTRIAL  RADIATION  30 

A.  THEORETICAL  AND  EMPIRICAL  BASIS  30 

B.  NET  FLUX  F10* 32 

1.  Net  Flux  F10*  with  Clear  Skies 33 

2.  Net  Flux  F10*  with  Overcast  Clouds  in 

Upper  Layer  Only  (k=4  to  k=6) 33 

3.  Net  Flux  F10*  with  Overcast  in  Lower 

Layer  Only  (k=8  to  k=9) 34 

4.  Net  Flux  F10*  with  Overcast  Clouds  in 

Both  Layers  (4,6)  and  (8,9) 34 

5.  Composite  F10*  Calculations  34 

C.  NET  FLUX  F6* 35 

D.  NET  FLUX  F2* 35 

E.  APPLICATIONS  TO  HEAT  BALANCE  COMPUTATIONS  —  36 

F.  COMPARISONS  WITH  PUBLISHED  RESULTS  37 

IV.  SOLAR  RADIATION 40 

A.  COMPOSITION  OF  SOLAR  INSOLATION  40 

B.  DISPOSITION  OF  F(S)  INSOLATION 42 

1.   Clear-Sky  Case 43 

5 


)ST 


2.  Cloudy-Sky  Cases  44 

3.  Composite  F(S)  Insolation 45 

C.  DISPOSITION  OF  F(A)  INSOLATION 45 

1.  The  Clear-Sky  Case 45 

2.  Overcast  in  Both  High  and  Low  Cloud 

Layers 47 

3.  Disposition  of  F(A)  Insolation  with  an 
Upper  Overcast 51 

4.  Disposition  of  F(A)  Insolation  with  a 

Low  Overcast 52 

5.  Weighted  F(A)  Layer-Absorptions  and 
Surface-Absorption  Insolation  55 

6.  Absorptivity  (ABA)  by  Layers  55 

D.  REFLECTION  OF  F(A)  AND  F(S)  ENERGY 57 

1.  Albedo 57 

2.  Composite  Absorptivity  (ABG)  by  the 
Earth's  Surface  58 

3.  Computational  Check  59 

E.  STATISTICAL  TESTS  59 

1.  Clear-Sky  Case  Regressions 59 

2.  Cloudy-Sky  Case  Regressions  61 

SENSIBLE  AND  LATENT  HEAT  TRANSPORT  AT  THE  SEA- 
AIR  INTERFACE 64 

A.  GENERAL  CONSIDERATIONS  64 

B.  EVAPORATION 64 

C.  SENSIBLE  HEAT  TRANSPORT 65 

D.  COMBINED  TURBULENT  HEAT  FLUX  ACROSS  THE 
SEA-AIR  INTERFACE  67 

6 


)NTER£Y, 


VI.  MERIDIONAL  CROSS-SECTIONAL  DEPICTION  OF  THE 
HEATING  MODEL  COMPUTATIONS  68 

A.  GENERAL  DESIGN  OF  HEAT-MODEL  OUTPUT  68 

B.  CROSS-SECTIONAL  DATA  AT  LEVEL  k=2 69 

C.  CROSS-SECTIONAL  RADIATIVE  TRANSFERS  IN 

LAYER  (2,6)  71 

D.  CROSS-SECTIONAL  DATA  FOR  LAYER  (6,10)  72 

1.  Radiative  Transfers  Only 72 

2.  Heat  Balance  at  the  Earth's  Surface  73 

3.  Sensible  Heat  Modification  of  the 
Atmospheric  Columns  73 

E.  MERIDIONAL  CROSS-SECTIONS  OF  THE  VERTICAL 

HEAT  BUDGET 74 

VII.  MEAN  BUDGETARY  BALANCE  MODEL  85 

A.  RADIATION  MODEL  OF  THE  OCEAN- ATMOSPHERE 

SYSTEM 85 

B.  RADIATION  BALANCES  OF  THE  OCEAN-TROPOSPHERE 
SYSTEM 89 

C.  TROPOSPHERIC  HEATING  92 

D.  SURFACE  HEATING  RATES  96 

VIII.  CONCLUSIONS 100 

LIST  OF  REFERENCES 102 

INITIAL  DISTRIBUTION  LIST  105 


NAVAL  POS 
MwNTHU.*, 


LIST  OF  TABLES 


Table 


I.  a.   Example  of  an  Original  FNWC  Gridpoint 

Sounding  for  Point  (1,1)  23 

b.   Example  of  the  Corresponding  Radiative 

Sounding  with  Temperature  and  Mixing  Ratio 

listed  at  k-levels 23 

II.  a.   Comparison  of  Net  Flux  at  the  Surface,  F10*, 

as  found  by  this  study  for  16  October  1973 

and  the  annual  mean  Fl 0  *  derivable  from 

Malkus  (1962,  Figs.  9,10)  39 

b.   Comparison  of  Net  Flux  to  space,  F2  ,  as  found 
by  this  study  and  by  Raschke  et  al .  (1973), 
based  on  NIMBUS  III  heat  budget  studies  39 

III.  Comparisons  of  Radiational  Balances  (a)  at  the 

tropopause  R  ,  (b)  in  the  vertical  column  R   be- 
s  a 

tween  k=2  and  k=10,  and  (c)  at  the  earth's  surface 
R  of  the  meridionally  averaged  values  as  found  in 
this  study,  the  annual  values  reported  by  Malkus 
(1962),  and  those  for  25  April  1973  found  by  Jenks 
(1974)  91 

IV.  Comparisons  of  Tropospheric  Radiative  cooling  rate 

(R  ),  turbulent  transport  warming  rate  (E+Hr),  and 
a  l 

net  tropospheric  heat  loss  rate  (Q   +S  )  as  found 
^  va   a 

by  this  study  for  16  October  1973  and  as  reported 

by  Malkus  (1962)  for  yearly  climatology  94 

V.  Comparisons  of  Net  Radiative  warming  rate  at  the 
surface  (R),  the  heat  loss  rate  at  the  surface  due 
to  turbulent  transport  -(E+Hr),  and  the  net  warm- 
ing or  cooling  at  the  earth's  surface  (Q   +S)  for 
16  October  1973  and  the  annual  mean  values  re- 
ported by  Malkus  (1962)  98 


8 


VAL  POS 
)NTEa£Y, 


LIST  OF  FIGURES 
Figure 

1.  Five-Layer  Radiative  Sounding  used  in  this 

study 18 

2.  FNWC  Polar  Stereographic  grid  and  meridians 

(lines  1,2,3,  and  4)  Selected  for  Study 21 

3.  Schematic  Representation  of  F(A)  Insolation 
Disposition  in  the  case  of  Two  Overcast  Layers  --   48 

4.  Schematic  Representation  of  F(A)  Insolation 
Disposition  with  an  Upper  Overcast  Layer  Only  53 

5.  Schematic  Representation  of  F(A)  Insolation 
Disposition  with  a  Lower  Overcast  Layer  Only  56 

6.  Key  to  Meridional  Cross-Sections  for  Figs.  7,..., 

10 76 

7.  125°W  Longitudinal  Cross-Section 

a.  Tropical  Section  77 

b.  Higher  Latitude  Section  78 

8.  170° W  Longitudinal  Cross-Section 

a.  Tropical  Section  79 

b.  Higher  Latitude  Section  80 

9.  145°E  Longitudinal  Cross  Section 

a.  Tropical  Section  81 

b.  Higher  Latitude  Section  82 

10.  35°W  Longitude  Cross-Section 

a.  Tropical  Section  83 

b.  Higher  Latitude  Section  84 


NAVAL  POS 
MONT  EOi.1, 


Figure 

11.  a.   Key  to  Mean  Radiation  Meridional  Cross- 

Section  for  Fig.  11(b)  and  (c)  86 

b.  Mean  Radiation  Meridional  Cross-Section  for 
Tropical  Latitudes  87 

c.  Mean  Radiation  Meridional  Cross-Section  for 
Mid-Latitudes  88 

12.  Radiational  Balance  at  the  Tropopause  (R  ),  in 
the  Atmospheric  column  between  k=2  and  k=10  (R  ), 

and  at  the  surface  (R)  for  16  October  1973  90 

13.  Tropospheric  Heat  Budget  Disposition  93 

14.  Surface  Heat  Budget  Disposition  97 


10 


NAVAL  POST 


LIST  OF  SYMBOLS 


A(m, n) 

a(m, n) 

a* 

ABA 

ABG 

ALB 

ATRAN 

BALB 


BALk^ 


BALT 
C 

CL 
CL(I) 

E 

e 
x 

F(A) 

FADJ 
Fkik2 


k 


solar  insolation  absorbed  in  the  layer  (m,n) 
Manabe-Moller  absorptivity  function 
turbulent  transfer  coefficient 
absorptivity  of  the  troposphere 
total  insolation  absorbed  by  earth's  surface 
earth-atmospheric  system  albedo 
transmissivity  of  the  troposphere 
Stefan-Boltzmann  blackbody  flux  at  T 
24-hour  averaged  radiational  balance  at  earth's 
surface 

24-hour  averaged  radiational  balance  for  layer 
(k!,k2) 

24-hour  averaged  radiational  balance  at  tropopause 
carbon  dioxide  layer  absorber  mass 
total  opaque  cloud  cover 

fractional  cloud  amount  for  layer:  1  =  1  in  600  to 
400  mb;  1  =  2  in  900  to  800  mb 
East  longitude;  evaporation 

vapor  pressure  at  top  of  constant  flux  layer 
solar  insolation  subject  to  water  vapor  absorp- 
tion only 

total  incoming  insolation  at  top  of  atmosphere 
net  infrared  flux  divergence  between  level  ki  and 
k2 
net  infrared  flux  at  level  k 


11 


EY. 


FNWC      -  Fleet  Numerical  Weather  Central 

F(S)  -  solar  insolation  subject  to  Rayleigh  scattering 
only 

g         -  gravity 

h         -  hour  angle 

H  -  height  of  homogeneous  atmosphere;  24-hour  aver- 
aged hour  angle 

H         -  sensible  heat  transport 

I         -  abscissa  grid  location 

IA10(m,n)  -  solar  insolation  absorbed  at  surface  with  cloud 
condition  (m,n) 

IS10(m,n)  -  solar  insolation  at  surface  subject  to  Rayleigh 
scatter  with  cloud  condition  (m,n) 

J         -  ordinate  grid  location 

k         -  pressure  level  used  in  this  study  equal  to  10a 

K*        -  eddy  turbulent  transfer  coefficient 

L         -  latent  heat  of  vaporization 

M         -  water-vapor  mass  path  length 

N         -  North  latitude 

P,         -  pressure  in  millibars  (mb)  at  level  k 

q,         -  mixing  ratio  at  level  k 

q  ,        -  saturated  mixing  ratio  at  level  k 
sk 

QAVE      -  24-hour  averaged  insolation  at  the  tropopause 

r         -  Bowen  ratio 

R         -  correlation  coefficient;  net  radiation  balance 

at  the  surface 

R         -  mean  radiative  cooling  rate  in  troposphere 
a 

12 


NAVAL  POST 
WwfiTEflti, 


Rd 

REF 
RE  FA 
REFS 
R.H. 
R_ 


k 

TRAN 


U 

W 

W(m,n) 

Z 

Zk 
a(G) 

a(R) 

6 
ewc 


A 

7T 


universal  gas  constant 

total  insolation  reflected  back  to  space 
F(A)  insolation  reflected  back  to  space 
F(S)  insolation  reflected  back  to  space 
relative  humidity 

mean  radiative  energy  gain  (loss)  rate  at  ocean- 
troposphere  system 

South  latitude;  effective  solar  constant;  heat 
storage  at  surface 

heat  storage  term  for  the  troposphere 
temperature  at  level  k 

total  insolation  incident  at  the  earth's  surface 
temperature  at  the  top  of  constant  flux  layer 
water-vapor  layer  absorber  mass 
West  latitude 

cloud  fractional  weight  for  cloud  condition  (m,n) 
Zenith  angle 
height  at  level  k 
surface  albedo 
Rayleigh  clear  sky  albedo 
critical  lapse  rate 
solar  declination  angle 

emissivity  due  to  water  and  carbon  dioxide 
absorber  mass  at  indicated  layer 
potential  temperature  at  level  k 
longitude 
surface  pressure;  pi  =  3.1416 


13 


5S 


p         -  density 

a         -  sigma  pressure  level  used  by  FNWC,  normalized 

to  surface  pressure 
d)         -  latitude 


14 


ACKNOWLEDGEMENT. 

The  author  wishes  to  express  his  appreciation  to  his 
thesis  advisor,  Professor  F.  L.  Martin,  for  his  suggestions, 
advice,  guidance  and  support  in  this  research. and  to  Mr. 
Russel  D.  Schwanz  for  his  expert  programming  assistance. 

Appreciation  is  also  expressed  to  the  author's  wife  for 
her  support  and  perseverance  and  to  Mrs.  Donna  Grigsby  for 
typing  this  manuscript. 


15 


I.   INTRODUCTION 

This  study  was  part  of  an  ongoing  project  to  produce  an 
empirical  heat  budget  parameterization  for  utilization  in  the 
Fleet  Numerical  Weather  Central  (FNWC)  primitive  equation 
prediction  model.   The  project  will  eventually  examine  the 
heat  balance  of  the  earth-ocean  system  using  representative 
mid-seasonal  data  from  each  of  the  four  seasons  of  the  year 
1973-74.   This  particular  study  covers  15-16  October  1973, 
using  gridded  FNWC  temperature  and  humidity  fields,  at  con- 
stant pressure  levels,  for  this  period. 

The  parameters  which  influence  the  radiative  calculations 
most  strongly  are  the  two  cloud  layers  of  fractional  amounts 
CL(1)  and  CL(2),  at  600  and  900  mb,  respectively,  at  each 
grid  point.   The  fractional  amounts  are  specified  using  the 
large-scale  formulations  developed  by  Smagorinsky  (1960)  as 
also  was  done  by  Jenks  (1974).   The  possibility  of  two  cloud 
layers  introduces  substantial  complexity  to  the  radiative 
physics  compared  to  the  one-layer  cloud  model  previously  de- 
scribed by  Martin  (1972).   The  adaptation  of  the  radiational 
model  to  FNWC  gridpoint  soundings  is  based  on  an  unpublished 
manuscript  (Martin,  1974).   The  Martin  manuscript  in  turn 
has  similarities  to  the  multi-cloud  heating  packages  of 
Arakawa  (1972)  and  Rodgers  (1967). 

Heat-balance  computations  were  made  at  the  ocean  surface, 
at  the  tropopause,  and  for  the  atmospheric  layer  surface  to 
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200  mb.   A  major  purpose  was  to  determine  whether  the  cloud- 
cover  parameterization  previously  adopted  by  Jenks  (1974) 
lends  itself  to  realistic  budgetary  computations  (as  compared 
to  the  results  of  others)  for  Oct.  15-16  over  the  Northern 
Hemispheric  oceans.   In  studying  the  net  heating  effects  at 
the  ocean  surface,  the  evaporative  and  sensible-heat  trans- 
fers were  modeled  after  Kesel  and  Winninghoff  (1972),  and 
Kaitala  (1974),  but  with  some  modifications  to  be  described 
later  (Section  V). 

To  assess  the  effects  of  the  cloud  cover  parameterization 
on  the  heat  transfers  of  this  study,  gridpoint  data  along 
four  meridians  were  selected  from  the  FNWC  analyses  of  15-16 
October  1973.   Three  of  the  meridians  are  from  the  Pacific 
Ocean  area  and  the  fourth  is  from  the  Atlantic  Ocean  (see 
Fig.  2).   Most  of  the  gridpoint  soundings  sampled  were  in 
the  Northern  Hemisphere. 

The  five  sigma  (a  =  p/tt)  analysis  levels  of  the  FNWC 
primitive  equation  model  were  adapted  for  data  input  in  the 
vertical  scale  of  this  study.   The  k-levels  (k=l , . . . , 10) , 
where  k=10a,  are  used  for  the  purpose  of  simplified  refer- 
encing of  the  data  levels  in  the  equations  derived  and  are 
shown  in  Fig.  1.   The  cloud  layers,  when  they  exist,  are  con- 
sidered to  be  in  two  layers,  a  high-level  cloud  layer  from 
k=4  to  k=6  and  a  low-level  cloud  deck  from  k=8  to  k=9,  with 
amounts  CL(1)  and  CL(2),  respectively.   In  pressure-scaling 
the  water  vapor  and  C02  absorber  masses  a  realistic  method, 
after.  Smith  (1966),  was  used  to  extrapolate  water  vapor 
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Figure  1.   Five-layer  radiative  sounding  used  in  this  study. 
Levels  are  identified  by  their  values  on  the  k-scale,  while 
layers  are  identified  by  their  level  boundary  indices  in 
parent heses ,  e.g.  (8,10).   Pressure-scaled  water  vapor  and 
C02  mass  increments  U  and  C,  respectively,  are  integrated 
with  respect  to  the  surface  and  introduced  at  even  levels 
while  the  mixing  ratio,  q,  is  formulated  at  odd  levels.   The 
temperature,  T,  is  expressed  for  all  levels. 
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mixing  ratio  values  above  the  level  of  k=4 ,  where  radiosonde 
moisture  data  is  normally  not  available.   The  net  flux  of 
terrestrial  radiation  computations  make  use  of  empirical  ab- 
sorptivities  of  water  vapor  and  C02  (Sasamori,  1968),  which 
are  incorporated  in  the  Martin  (1974)  manuscript.   The  solar 
radiation  package  in  Martin  (1974)  has  been  modeled  after 
Arakawa  (1972),  except  that  only  two  cloud  layers  are  con- 
sidered possible,  as  contrasted  with  three  in  the  Arakawa 
version. 
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II.   DATA  PREPARATION 

The  data  used  in  this  study  was  in  the  form  of  gridpoint 
soundings  chosen  along  four  meridians  of  FNWC  63-by-63  grid- 
mesh  temperature  and  humidity  data  fields  to  insure  minimum 
over  land  data.   The  three  meridians  (and  their  respective 
number  of  soundings)  selected  over  the  Pacific  Ocean  were 
located  at  125°W  -  (25  soundings),  170°W  -  (25  soundings), 
and  145°E  -  (17  soundings).   The  Atlantic  Ocean  meridian  was 
35°W  -  (26  soundings).   Fig.  2  shows  how  the  meridians  were 
located  on  the  FNWC  polar  stereographic  grid.   This  procedure 
of  choosing  data  lines  along  diagonals  or  along  J=32  (through 
the  pole-point)  of  the  FNWC  data-fields  obviates  the  neces- 
sity of  interpolating  spatially  to  gridpoints.   Data  along 
line  3  in  the  Pacific  were  not  extended  southward  of  (9,55) 
because  they  fell  over  land  masses  (New  Guinea  and  Northern 
Australia)  where  the  surface  temperatures  were  too  high  to 
be  representative  of  the  oceanic  values. 

The  vertical  soundings  on  these  data  lines  were  taken 
from  the  original  FNWC  63-by-63  gridpoint  analysis  of  tem- 
perature at  the  surface  and  at  nine  other  standard  levels  up 
to  100  millibars  for  0000GMT,  16  October  1973  for  the  Pacific 
Ocean  data  cases  and  12  hours  earlier  for  the  Atlantic  Ocean 
data  cases.   In  addition  to  these  standard  level  values,  a 
T   value  was  listed  which  is  the  FNWC  calculated  temperature 
(see  Kesel  and  Winninghoff,  1972)  for  the  top  of  the  turbu- 
lent layer  of  constant  sensible  heat  flux.   This  level, 
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Figure  2.   FNWC  polar  stereographic  grid  and  meridians  (lines 
1,  2,  3,  and  4)  selected  for  study.   The  longitudes  A   are 
shown  for  each  meridian  as  well  as  the  extent  considered  of 
each  meridian. 
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which  is  approximately  20  meters  above  the  surface,  was  des- 
ignated by  a  code  label  of  9999  under  the  column  labeled 
pressure  (Table  1(a)).   The  six  moisture  values  of  this  table 
were  given  by  the  vapor  pressure  at  sea  level,  and  by  dew- 
point  depression  for  levels  from  925  mb  to  400  mb.   In  addi- 
tion, the  parameter  e  was  provided  by  FNWC  as  the  vapor 
pressure  at  the  top  of  the  constant  flux  layer  at  the  level 
indicated  by  the  code  9999  (Table  1(a)).   Table  1(a)  gives 
an  example  of  the  original  sounding  for  gridpoint  (1,1)  of 
the  FNWC  map  for  the  given  period. 

In  order  to  perform  radiative  calculations,  it  is  neces- 
sary to  have  water  vapor  and  C02  absorber  masses  and  cloud 
amounts  within  the  required  k-level  boundaries  (Fig.  1).   To 
do  this,  a  systematic  method  was  applied  to  arrive  at  the 
temperature  and  mixing  ratio  distributions  (Table  1(b))  as 
detailed  in  the  following  subsections.   It  should  be  noted 
that  all  soundings  start  at  sea  level,  but  as  a  reasonable 
approximation  surface  pressure  was  considered  to  be  tt  =  1000  mb 
Therefore  the  k-levels  of  Fig.  1  become  1000.,  900.,  800., 
....,  200.,  100.,  0.0  mb,  respectively. 
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Table  1(a).   Example  of  an  original  FNWC  gridpoint  sounding 
for  point  (1,1).   The  humidity  parameters  for  the  surface 
and  the  tope  of  the  constant  flux  layer  (9999)  are  vapor  pres- 
sures and  those  between  925,..., 400  mb  are  dew-point  depres- 
sion. 


Pressure  (mb) 

surface 
925 
850 
700 
500 
400 
300 
250 
200 
150 
100 
9999* 


T(°C) 


Humidity  Parameters 


25.5 

25.6  mb 

19.3 

3.6°C 

15.5 

3.2°C 

6.5 

5.4°C 

-10.2 

5.4°C 

-21.9 

2.1°C 

-36.8 

-45.8 

-56.0 

-67.3 

-80.6 

T  =  24.7 

X 

e   =  26.3  mb 

X 

Table  1(b).   Example  of  the  corresponding  radiative  sounding 


with  temperat 

,ure 

and 

mixing  rat 

io 

listed 

at 

k-levels . 

Pressure 

(mb) 

T(°K) 

M: 

Lxing  Ratio  (g/kg) 

1000 

298.7 

15.92 

900 

290.9 

11.45 

800 

285.9 

700 

279.7 

5.95 

600 

272.4 

500 

263.0 

2.31 

400 

251.3 

1.42 

300 

236.4 

0.56 

200 

217.2 

100 

192.6 

0.03 

0 

192.6 

0.00 

*  Values  at  top  of  constant  flux  layer  as  used  in  FNWC  heat- 
ing package  (Kesel  and  Winninghoff,  1972). 
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A.  TEMPERATURE  TERMS 

The  temperatures  were  listed  at  each  mandatory  level  be- 
tween the  surface  and  100  mb .   The  temperature  Tl0    was  set 
equal  to  the  listed  surface  air  temperature.   At  the  top  of 
the  atmosphere,  the  temperature  was  assumed  to  be  isothermal 
from  100  mb  to  P=0  mb,  with  the  mixing  ratio  approaching 
zero  at  P=0  mb.   The  radiative  sounding  temperatures  for  the 
remaining  k-levels  were  obtained  from  either  their  corre- 
sponding listed  temperature  or  by  a  three-point  Lagrangian 
interpolation  scheme  of  the  form: 

Tk  =  To((Pk-Pi)(Pk-P2))/((Po-Pi)(Po-P!)) 

+  T1((Pk-P0)(Pk-P2))/((P1-P0)(Pi-P2))  (2-1) 

+  T2((Pk-P0)(Pk-P1))/((P2-Po)(P2-Pi)) 

Here  k  is  the  pressure  level  desired  for  T,  interpolation, 
P0  lies  above  P,  ,  and  Px ,  P2  lie  below  level  k.   The  param- 
eters  T0 ,  Tx ,  and  T2  are  temperatures  at  the  same  levels  as 
P0 ,  Pi  and  P2 ,  respectively.   As  an  example,  if  k=9,  Pi  and 
P2  were  taken  as  925  mb  and  1000  mb,  respectively. 

The  T   temperature  was  not  utilized  in  the  interpolation 
scheme  described  above,  but  was  carried  along  as  a  twelfth 
temperature  for  later  use  in  the  heat  budget  computations. 

B.  MOISTURE  TREATMENT  IN  THE  SOUNDING 

The  moisture  parameters  from  the  original  FNWC  soundings 

(Table  1(a))  were  converted  into  mixing  ratios  at  each  of 

the  original  sounding  levels.   At  k=10  the  surface  vapor 

pressure,  e  .  ,  was  used  to  calculate  the  mixing  ratio  as 
air 
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follows: 

q100o  *  .62197  e    /1000  (2-2) 

all 

where  e  .   is  in  millibars.   To  obtain  a  mixing  ratio  for 
air 

the  top  of  the  constant  flux  layer,  Equation  (2-3)  was  used 
to  compute  q   from  e 

A.  X 

q   b  (e  /1000)  .62197  (2-3) 

-A.  X 

where  e   is  the  vapor  pressure  at  the  top  of  the  constant 
flux  layer  in  mb. 

The  mixing  ratio  at  level  P  is  calculated  at  the  remain- 
ing original  data  levels  by  the  use  of  Equation  (2-4)  which 
is  taken  from  Fleagle  and  Businger  (p.  63,  1963); 

621  97  e(A-B/T)      "B(T-T  ) 
q(P)  =  621.97  e exp[       D  ]  (2-4) 

UD 
A   =  21.656 
B   =  5418. 0°K 

T   =  temperature  °K  at  level  P 
T   =  dew  point  (°K)  at  level  P 

The  calculated  q-values  were  then  interpolated  to  the  k-lev- 
els  using  an  equation  similar  to  Equation  (2-1).  An  example 
of  the  interpolation  results  for  gridpoint  (1,1)  is  shown  in 
Table  1(b). 

Since  accurate  humidity  information  is  not  generally 
available  for  the  upper  troposphere,  a  logarithmic  extrapola- 
tion formula  was  used  to  find  q-values  at  k=l,  2,  and  3. 

aIZi  =  (_^-)A  (2-5) 

q5     ^500; 
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This  formula  as  described  by  Jenks  (1974)  is  similar  to  that 

developed  by  Smith  (1966).   It  is  solved,  using  a  "best-fit" 

a  P 

between  variables  y  =  log  -*—  and  x  e  log  tttt:,  for  the  param- 

q  5  OUU 

eter  A.   The  six  q-values  of  the  original  soundings  corre- 
sponding to  P  =  1000,  925,  850,  700,  500,  and  400  mb  are 
used  to  determine  X,  which  is  then  used  in  Equation  (2-5)  to 
extrapolate  the  q  values  to  P  =  300,  200,  and  100  mb. 

A  correlation  coefficient  was  statistically  computed  in 
connection  with  each  A-profile  of  the  form: 


6 

£  x2 

The  correlation  coefficient  enabled  the  monitoring  of  the 
usefulness  of  the  A-profile  technique  and  was  generally 
found  to  be  in  the  range  of  0.95  to  0.99. 

C.   PRESSURE-SCALED  ABSORBER  MASSES 

Having  computed  the  mixing  ratios  at  each  odd  k-level, 
it  is  now  possible  to  compute  the  pressure-scaled  water  va- 
por absorber  masses  in  each  layer.   In  Figure  1,  these  water 
vapor  layer-masses  were  denoted  by  U,  however  here  we  will 
use  notation  which  is  more  representative  of  the  particular 
layer  under  consideration.   For  example,  the  pressure-scaled 
water  vapor  mass  in  the  lowest  layer  (8,10)  is  given  by 

UC8.10)  -3|^  t^-^]0'72  (2-6) 


26 


where  the  Moller  and  Raschke  (1964)  pressure-scaling  exponent 
0.72  has  been  employed.   For  the  layer  (6,8),  an  analog  to 
Equation  (2-6)  was  used  by  replacing  the  parameters  with  the 
subscript  9  with  their  counterpart  parameters  with  the  sub- 
script 7.   A  similar  procedure  was  followed  for  the  remain- 
ing layers  (4,6),  (2,4),  (0,2).   In  each  of  the  five  layers 
AP  =  200  mb. 

The  algorithm  for  computing  the  integrated  water  vapor 
mass  above  the  earth's  surface  is  as  follows: 

U(10)  =  0.0 

U(8)  =  U(10)  +  U(8,10) 

U(6)  =  U(8)  +  U(6,8) 

U(4)  =  U(6)  +  U(4,6)  (2-7) 

U(2)  =  U(4)  +  U(2,4) 

U(l)  =  U(2)  +  0.5U(0,2) 

U(0)  =  U(2)  +  U(0,2) 

A  similar  computation  scheme  is  used  for  the  carbon  di- 
oxide scaled  mass.   Here,  however  it  is  customary  to  repre- 
sent the  pressure-scaled  mass  in  terms  of  the  normal 
temperature  and  pressure  (N.T.P.)  volume  over  a  cm2  of  the 
earth's  surface.   Thus  for  the  lower  layer: 

C(8,10)  =  3.14  x  10--(|g-)(101^25)°'72  (2-8) 

°  s 

The  fractional  amount  3.14  *  10-1*  is  the  constant  mixing  ra- 

AP 
tio  of  C02  by  volume,  and  is  the  N.T.P.  thickness  of  the 

gPs 

layer  AP  =  P10-Pe.   The  quantity  p   is  constant  density  in  a 
standard  homogeneous  atmosphere  of  depth  H  where  T  =273.16°K 

27 


ai 


everywhere  and  P   =  1013.25  mb.   The  parameter  H,  also  gen- 
erally  known  as  the  scale  height  of  the  homogeneous  atmos- 
phere, may  be  further  defined  as: 

RT 

H  =  — -  =  7.995  x  io5  cm 
g 

Therefore  Equation  (2-8)  may  be  rewritten  as 

AP  P       0.72 

C(8,10)  =  3.14  x  io-  HC-ff-^Kj^gg)         (2-9) 

All  other  layer  thicknesses  of  C02  are  obtained  in  a  similar 
manner  to  Equation  (2-9)  by  substituting  the  appropriate 
subscripted  indices  7,5,...,1  on  the  pressure  P. 

The  algorithm  for  computing  the  integrated  carbon  diox- 
ide mass  above  the  earth's  surface  is  as  follows: 

C(10)  =  0.0 

C(8)   =  C(10)  +  C(8,10) 

C(6)   =  C(8)  +  C(6,8) 

C(4)   =  C(6)  +  C(4,6)  (2-10) 

C(2)   =  C(4)  +  C(2,4) 

C(l)   =  C(2)  +  0.5C(0,2) 

C(0)   =  C(2)  +  C(0,2) 

D.   CLOUD  PARAMETERIZATION 

In  order  to  determine  the  fractional  cloud  cover  amounts 
in  the  layers  (4,6)  and  (8,9),  the  saturated  vapor  pressure 
must  be  calculated  at  levels  k=5  and  k=9,  respectively  by 
the  following  equations: 

fi?1Q7  P(A-B/T5> 
qs(500)  =  '621975^Q (2-lla) 
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62197  e(A~B/T9) 
qs(900)  =  -621979^0 (2-llb) 

where  A  =  21.656  and  B  =  5418. 0°K  as  before  in  Eq.  (2-4). 

Knowing  the  saturated  vapor  pressure  at  levels  k=5  and 
k=9  enables  the  determination  of  the  relative  humidity  at 
those  levels: 

R-H-<9>  =ftf§M)  (2"12b) 

The  relative  humidities  were  then  substituted  into 
Smagorinsky ' s  (1960)  formulation  for  fractional  cloud  amounts 

CL(1)  =  2.0  (RH(5))  -  0.7  (2-13a) 

CL(2)  =  3.33  (RH(9))  -  2.0  (2-13b) 

These  formulations  were  employed  by  Jenks  (1974)  and  are 
also  to  be  tested  in  this  study  in  which  a  different  season- 
al sample  is  considered. 

These  fractional  cloud  amounts  are  a  function  of  large 
scale  influences  only  and  do  not  consider  small  scale  con- 
vective  activity,  seasonal  conditions  or  latitude.   The  use 
of  such  clouds  is  for  radiational  study  only. 
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III.   TERRESTRIAL  RADIATION 

A.   THEORETICAL  AND  EMPIRICAL  BASIS 

Empirical  formulas  were  developed  by  Sasamori  (1968)  for 
computing  emissivities  for  long-wave  flux  calculations  asso- 
ciated with  the  NCAR  General  Circulation  Model.   These  em- 
pirical emissivity  formulas  can  be  related  to  the  values 
built  into  Yamamoto's  (1952)  well  known  Radiation  Chart, 
which  was  developed  from  radiative  transfer  theory.   This 
chart  has  proved  to  be  very  accurate  for  computational  pur- 
poses, and  was  used  by  Jenks  (1974)  as  a  schematic  guide  for 
integration  of  the  radiative  transfers  through  the  various 
layers  of  interest  in  the  soundings  discussed  in  Section  II. 

The  necessary  emissivity  formulas  of  Sasamori  for  utili- 
zation in  the  FNWC  heating  package  are  those  which  will  yield 
expressions  for: 

F10*  =  Net  IR  flux  at  earth  surface 

F6*   =  Net  IR  flux  at  level  k=6 

F2*   =  Net  IR  flux  at  level  k=2 

F610  =  Net  IR  flux  divergence  in  the  layer  (6,10) 

F26   =  Net  IR  flux  divergence  in  the  layer  (2,6) 

A  workable  scheme  for  making  these  computations  is  presented 
below  for  various  combinations  of  cloud  cover  [CL(1),  CL(2)]. 
For  example,  in  order  to  compute 

Bio 

F10*  -  Bio  -  B£0  ewc(U,C,T)dB  (3-1) 
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where  Bk  is  the  Stef an-Boltzmann  blackbody  flux  at  T  ,  one 

must  have  a  good  representation  of  the  emissivity  e    as  a 

wc 

function  of  both  water  vapor  and  carbon  dioxide  absorber 

masses  (U  and  C,  respectively)  along  the  sounding  which  is 

the  path  of  integration.   For  the  integral  evaluation  of 

Eq.  (3-1),  Sasamori  (1968)  proposed  temperature-independent 

emissivity  formulations  z        as  exemplified  for  the  absorber 

wc        ^ 

masses  in  layer  (8,10)  as  follows: 

e   (8,10)  =  {.240  log10[U(8,10)]+.622}+.07262{(l. 

WO 

-. 62556 [U( 8 ,10) +.0286]' 26 ) [ logi 0C(8, 10) 

+1.064]}  (3-2) 

Formulas  of  the  type  given  by  Eq .  (3-2)  apply  for  T  >  210° K. 
The  quantity  in  the  first  brace  of  Eq.  (3-2)  indicates  the 
water  vapor  emissivity  alone  and  that  in  the  second  brace  is 
the  additional  contribution  by  that  part  of  the  C02  emissiv- 
ity not  already  included  in  the  water  vapor  emissivity  ex- 
pression . 

Eq.  (3-2)  was  used  in  the  calculation  of  the  atmospheric 
emissivities  over  all  layers  (6,10),  (4 , 10) , . . . . , (1 , 10) ,  for 
which  cumulative  mass  values  U(6,10),  C(6,10),  etc.,  were 
available  from  the  radiative  soundings.   Thus  Eq.  (3-2) 
proved  very  versatile.   For  completion  of  the  integration 
from  level  k=l  to  the  origin  (B=0,  e=1.0)  of  Yamamoto's  Ra- 
diation Chart,  the  U(T  ) ,  C(T  )  soundings  were  taken  to  be 

K  K 

isothermal  at  T=Tj  to  the  top  of  atmosphere  (to  level  "0"). 
This  leaves  only  the  temperature-dependent  emissivity  e 
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[U(0,10),T1]  which  applied  for  T  <  210°K,  also  computed  after 
Sasamori  (1968)  by  means  of 

I        UCCIO),^   =  8.34Tl[-353  log10U(0,10)-.44] 

wc  * 

x  U(0,10)[~'03455  logioU(0,10)-.705] 

x  [8.0/(.353  log10U(0,10)+3.56)] 

+  .07262{ ( 1. -. 62556 [U( 0,10) +.0286] ' 26 ) 

x  [logi0C(0,10)+1.064]}  (3-3) 

The  first  term  is  temperature-dependent,  as  well  as  dependent 
upon  the  summed  absorber  masses  relative  to  the  reference 
level  (in  this  example,  level  k=10).   The  term  in  the  brace 
involves  essentially  the  same  form  of  the  C02  correction  to 
the  total  emissivity  as  seen  in  Eq.  (3-2)  except  that  it  re- 
fers to  the  full  layer  (0,10).   Eq.  (3-3)  is  applicable  with 
appropriate  final  parameters  U(0,k),  C(0,k),  and  Ti  to  the 
determination  of  the  downward  flux  contributions  through 
level  k  for  all  of  the  reference  level  cases  to  be  consid- 
ered below,  i.e.  k=2,6,10. 

B .   NET  FLUX  F i „ * 

The  radiative  soundings  were  evaluated  in  the  form  of 
the  parameters  U(k,10),  C(k,10),  and  T,  ,  level  by  level. 

K 

Also  listed  for  each  radiative  sounding  are  the  cloud  cover 
parameters  CL(1)  and  CL(2),  which  in  general  are  non-zero. 
The  grid  area  may  be  visualized  as  composed  of  areal  frac- 
tions or  weights  defined  as  follows: 
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W(0,0)  =  (1-CL(1))(1-CL(2))  :  totally  clear  fraction 

W(1,0)  =  (CL(1))(1-CL(2))  :  overcast  in  upper  layer  on]y 

(3-4) 
W(0,1)  =  (1-CL(1))(CL(2))  :  overcast  in  lower  layer  only 

W(l,l)  =  (CL(1))(CL(2))  :  overcast  in  both  layers 

The  composite  radiative  effect  of  the  true  cloud  cover 
over  the  grid  area  is  then  reconstructed  by  multiplying  the 
results  of  the  respective  cases  defined  above  by  the  appro- 
priate weight  factors.   The  four  cases  will  now  be  defined 
in  order. 

1.   Net  Flux  F1Q*  with  Clear  Skies 

Defining  the  Stef an-Boltzmann  blackbody  flux  at  T 
as  B  ,  where  B   =  1.170403  x  10~7  Tj*  ly(day)-1,  the  net  flux 
at  the  surface  under  the  clear  sky  case  is  defined  as  fol- 
lows : 

F10*(O,O)  =  B10-.5{ewc(8,10)(B10-B8)+(B8-B6)[ewc(6,10) 


(3-5) 


+e   (8,10)]+[e   (4,10)+e   (6 , 10) ] (B6-B4 ) 
wc      /J  L  wc        wc 

+(B,-B2)[ewc(2,10)+ewc(4,10)]+[ewc(l,10) 
+ewc(2,10)](B2-B1)+ewc[(0,10),T1]B1} 

2.   Net  Flux  F10*  with  Overcast  Clouds  in  Upper  Layer 
Only  (k=4  to  k=6) 


F10*(1,0)  =  (B10-Be){l.-.5[ewc(8,lO)(B10-B8) 
+(ewc(8,10)+ewc(6,10))(B8-B6)]} 


(3-6) 
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3.  Net  Flux  F1Q*  with  Overcast  in  Lov/or  Layer  Only 
(k=8  to  k=9) 

F10*(0,1)  =  (B10-B9)[l.-.5ewc(9,lO)]  (3-7) 

4.  Net  Flux  FIQ*  with  Overcast  Clouds  in  Both  Layers 
(4,6)  and  (8,9) 

F10*(1,D  =  (B10-B9)[l.-.5ewc(9,10)]  (3-8) 

It  should  be  noted  that  Eq.  (3-7)  and  Eq.  (3-8)  are  identi- 
cal. 

5.  Composite  F10*  Calculations 

The  composite  results  are  obtained  by  multiplying 
the  above  defined  four  cases  by  their  respective  weights  as 
defined  by  Eq.  (3-4)  to  obtain 

F10*  =  [l-CL(2)]{(B10-B6)-.5cwc(8,10)(B10-B8) 

+[ewc(8,10)+ewc(6,10)](B8-Be)} 

+(l-CL(2))(l-CL(l)){B6-.5[(ewc(6,10) 

+ewc(4,10))(B6-Blt)  +  (ewc(4,10)+ewc(2,10))(B,-B2) 

+(ewc(2,10)+ewc(l,0))(B2-B1)+ewc((0,10),T1)B1]} 

+CL(2){(B10-B9)[l-.5ewc(9,lO)]}  (3-9) 

It  is  clear  that  Eq .  (3-9)  reduces  to  F10*(O,O)  when  CL(1)  = 
CL(2)  =  0,  and  to  F10*(l,l)  when  CL(1)  =  CL(2)  =  1.0.   In  a 
similar  manner  the  other  two  cloud  cases  may  be  recovered  by 
setting  the  weights  to  CL(1)  =  1.0,  CL(2)  =  0  and  CL(2)  =  1.0, 
CL(1)  =  0,  respectively. 
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The  above  described  procedure  for  obtaining  the  net 
flux  can  similarly  be  extended  to  levels  k=6  and  k=2  as 
listed  in  the  following  subsections. 

The  composite  Fi0*  values  will  be  compared  in  Section 
III.F  with  climatology  (Malkus,  1962)  to  ensure  the  reason- 
ableness of  this  calculation. 

C.  NET  FLUX  F6* 

The  individual  cases  for  F6*  can  be  obtained  from  the 
composite  F6*  as  described  above  for  F10*.  The  numerical 
from  using  the  Sasamori  type  emissivities  [Eqs.  (3-2)  and 
(3-3)]  for  the  composite  F6*  is  as  follows: 

F6*  =  [l-CL(l)]{B8-.5[e   (6,8)(B8-B6)+e   (4 , 6) (B6-Bu ) 

+(ewc(4,6)+ewc(2,6))(B„-B2)+(ewc(2,6)+£wc(l,6))(B2-B1) 

+ewc((0,6),T1)B1]}+(l-CL(l))(l-CL(2)){(B10-B8)[l. 

-.5(ewc(6,8)+ewc(6,10))]}+CL(l){(B8-B6)[l-.5(ewc(6,8))]} 

+CL(l)(l-CL(2)){(B10-B8)[l-.5(e   (6,8)  +  e:   (6,10))]} 

WC        WC         (3-10) 

D.  NET  FLUX  F2* 

In  a  similar  fashion  to  F10*  and  F6*  the  net  flux  at 
level  k=2  can  be  calculated  from  the  radiation  chart  using 
the  trapezoidal  rule  as  described  by  Jenks  (1974).   Again 
the  numerical  formulation  of  the  individual  cloud  cases  can 
be  obtained  from  the  composite  net  flux  formula  as  described 
in  Subsection  III.B.5.   The  composite  F2*  formula  is 
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F2*  =  [l-CL(l)]{B8-.5[ewc(2,4)(B^B2)+(ewc(2,4) 

+ewc(2,6))(B6-BJ  +  (E:wc(2)6)+cwc(2)8))(B8-B$) 

+ewc(l,2)(B2-B1)+cwc((0,2),T1 )B1])+(1-CL(1))(1-CL(2)) 

x{(B10-B8)[l-.5(ewc(2,8)+Gwc(2,10)]}+CL(l){Bl4 

-.5[ewc(2,4)(BIt-B2)+ewc(l,2)(B2-B1)+ewc((0,2),T1)B1]} 

If  the  terms  e   ((0,2),T1)B1  and  e   (1,2)(B2-B1)  are  set 
to  zero  in  Eq.  (3-11)  a  new  parameter  can  be  defined  as  the 
net  outgoing  long-wave  radiation  and  labeled  F2 .   This  pa- 
rameter was  computed  to  compare  with  satellite  climatolog- 
ical  observations  over  the  oceans  (Raschke,  et  al.,  1973). 


E.   APPLICATIONS  TO  HEAT  BALANCE  COMPUTATIONS 

For  computation  of  the  radiative  "balance"  at  the  top  of 
the  earth-tropospheric  system  only  F2*  of  Eq.  (3-11)  and  the 
total  insolation  absorbed  below  k=2  are  required.   At  the 
surface,  a  radiative  balance  is  computed  at  each  gridpoint 
using  the  absorbed  insolation  and  the  net  flux  F10*  of  Eq . 
(3-9).   Also  at  the  surface  a  "convective"  balance  is  com- 
puted by  including  the  contributions  of  evaporation  and  sen- 
sible heat  transport  (see  Sec.  V).   Net  flux  F10*,  evaporation 
and  sensible  heat  transport  are  all  considered  negative  con- 
tributions at  the  surface.   The  atmospheric  radiative  bal- 
ance requires  knowledge  of  the  long-wave  cooling  effects 
caused  by  the  flux  divergences  F26  and  F610  as  defined  by: 

F26  =  F2*  -  F6*  (3-12a) 

F610  =  F6*  -  Fio*  (3-12b) 
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As  will  be  seen  in  a  later  section,  these  parameters  act  as 
destabilizing  influences  in  their  respective  layers  when  as- 
sessed in  relationship  to  other  processes  in  cross  sectional 
display  form  in  Section  VI. 

F.   COMPARISONS  WITH  PUBLISHED  RESULTS 

Table  II  gives  a  visual  comparison  of  the  outgoing  long- 
wave radiation,  F2  (as  contrasted  with  F2*),  found  by  this 
study  as  compared  to  the  mean  values  published  by  Raschke 
et  al.  (1973).   Also  displayed  are  this  study's  global  net 
flux,  F10*,  computed  using  Eq.  (3-9)  as  compared  to  the  an- 
nual mean  deduced  from  Malkus  (1962,  Figs.  9,10).   The  16 
October  1973  results  have  been  interpolated,  gridpoint  by 
gridpoint,  to  5°  latitude  increments  for  each  of  the  four 
lines  (see  Fig.  2)  and  then  latitudinally  averaged  to  arrive 
at  the  "global"  value  for  each  latitude.   The  units  of  the 
values  from  Malkus  (actually  due  to  Budyko  (1956))  have  been 
converted  from  kcal  cm-2yr-1  to  those  of  cal  cm-2min_1  used 
in  this  study.   The  global  mean  values  listed  at  the  bottom 
of  each  column  were  derived  by  cosine  weighting  the  latitu- 
dinal values  listed  in  Table  II  (see  also  Sec.  VII). 

It  should  be  noted  that  the  values  of  F2  and  F10*  found 
by  this  study  for  latitudes  20°S,  10°S,  and  60°N  are  averages 
of  only  two  "sounding"  lines  while  the  other  latitude  values 
are  averaged  over  all  four  meridional  lines,  the  difference 
being  due  to  the  number  of  soundings  available  near  the  ex- 
tremities of  the  FNWC  cross-sections. 
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The  consistently  smaller  values  of  net  flux,  Fi0*,  found 
by  this  study  might  be  explained  by  the  fact  that  Budyko's 
results  are  yearly  climatological  values,  whereas  the  Fio* 
values  of  this  study,  involving  data  from  16  October  1973, 
occurs  near  the  time  of  the  oceanic  maximum  of  temperature 
and  higher  than  annual  mean  cloudiness.   This  would  reduce 
the  net  flux  Fio*  computed  here.   A  latitudinal  maximum  net 
flux,  Fio*,  was  observed  around  40°N  both  for  this  study  and 
the  comparison  figures  of  Malkus.   At  this  location,  Figs. 
7(b),  9(b),  10(b)  show  that  there  was  a  minimum  opaque  cloud 
cover  for  the  time  of  this  study. 

The  comparison  of  outgoing  long-wave  fluxes,  F2  (Table 
11(b)),  appears  very  consistent  when  comparing  only  northern 
hemispheric  latitudes,  considering  that  Raschke's  values  are 
based  upon  the  mean  for  a  month  rather  than  that  for  a  single 
day.   Here  maximum  values  of  very  nearly  the  same  order  oc- 
cur in  the  subtropics  where  minimum  cloudiness  occurs  to  im- 
pede the  emission  of  outgoing  terrestrial  flux  from  the  earth 
and  lower  atmosphere.   There  is  in  both  this  study  and  in 
Malkus',  a  well  defined  poleward  gradient  of  F2  by  latitudes. 

The  comparison  of  the  computations  of  Fio*  and  F2  made 
here  with  climatology  (Table  Ila.b)  cannot  confirm  or  deny 
the  specific  cloud-parameterization  formulas  (Eq.  2-13a,b) 
used  in  this  study,  but  they  indicate  that  they  have  a  gen- 
eral reliability. 
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(a) 


LAT 

Fi  o  * 

Fi  o  * 

(Malkus) 

(This  Study) 

20°S 

.100 

.068 

10°  S 

.090 

.054 

0° 

.080 

.068 

10°N 

.090 

.079 

20°N 

.114 

.078 

30°N 

.110 

.081 

40°N 

.120 

.107 

50°N 

.110 

.086 

60°N 

.120 

.067 

Global 

Mean 

.1018 

.0757 

(b) 


F2 


(This  Study)   (Raschke  et  al.) 


.333 

.341 
.355 
.378 
.389 
.388 
.376 
.346 
.315 

.3601 


.40 
.39 
.36 
.34 
.40 
.37 
.34 
.31 
.29 

.3615 


Table  11(a).   Comparison  of  net  flux  at  the  surface,  Fi0*, 
as  found  by  this  study  for  16  October  1973  and  the  annual 
mean  F10*  derivable  from  Malkus  (1962,  Figs.  9,10). 

(b).   Comparison  of  net  flux  to  space,  F2 ,  as  found 
by  this  study  and  by  Raschke  et  al.  (1973),  based  upon 
NIMBUS  III  heat  budget  studies. 
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IV.   SOLAR  RADIATION 

A.   COMPOSITION  OF  SOLAR  INSOLATION 

At  the  top  of  the  atmosphere  (k=0)  this  study  assumed  a 
solar  constant  of  2.00  ly  min-1  (Joseph,  1971);  furthermore, 
this  constant  was  assumed  subject  to  a  four  percent  attenua- 
tion above  the  tropopause  due  to  ozone  and  oxygen.   Thus  the 
effective  solar  constant  at  level  k=2  in  this  study  is  1.92 
ly  min-1 . 

To  compute  the  effective  solar  insolation  at  the  tropo- 
pause the  following  formula  was  used 

F(2)  =  S[^-]~2  Cos  Z  (4-1) 

m 

where   S     =  effective  solar  constant  at  level  k=2 

Cos  Z  =  cosine  of  the  zenith  angle,  a  function  of 

the  Julian  date 

r/r   =  ratio  of  the  actual  earth-sun  distance  to 
'  m 

the  mean  earth-sun  distance,  a  function  of 
the  Julian  date. 

The  Smithsonian  Meteorological  Tables  (List,  1958)  list  the 

ratio  r/r   for  15  October  at  0000GMT  as  0.997165,  this  value 
'  m 

was  utilized  for  all  four  data  lines.   The  cosine  of  the 
zenith  angle  was  determined  by 

Cos  Z  =  Sin  <J>  Sin  6  +  Cos  <|>  Cos  6  Cos  h        (4-2) 
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where 

<f>  =  latitude 

6  =  solar  declination  for  15  October  =  -8°  13.5' 
(Table  169,  List,  1958) 

h  =  hour  angle 

The  value  of  sine  <j>  was  calculated  using  one  of  two  differ- 
ent formulae,  depending  on  the  data-line  used  for  the  compu- 
tations, in  terms  of  the  FNWC  map  coordinates  (I, J)  as  in 
Eq.  (4-3a,b),  conversely  for  these  lines  one  may  solve  for  I 
in  terms  of  sine  $    as  in  (4-3c,d): 

Lines  1,3,4    Sin  <p    =    jffij;^  ;  2(32-1)*  (4"3a) 

t-  o  q-       a         973.752    -    (32-I)2  fA    Q.* 

Line   2  Sin  *   =   973>752   +    \32-iy  (4"3b) 

Lines  1,3,4    I  =  32  -  22.065  [■=  ^.^  .]  (4-3c) 

'  '  1  +  Sin  <j> J 

Line  2         I  =  32  -  31.205  [±   +°gin  ,]  (4-3d) 

for  Lines  1,2 
for  Line  3 
for  Line  4 

Here  I  is  the  abscissa  distance  on  the  FNWC  grid,  Fig.  2, 
and  varies  by  line  as  described  in  Sec.  II.   The  soundings 
for  lines  1,  2,  and  3  were  all  taken  at  0000GMT  with  the  so- 
lar  noon  existing  at  the  180   meridian,  therefore  the  hour 
angles  for  these  three  lines  were  55° ,  10° ,  and  35° ,  re- 
spectively.  For  line  4,  the  soundings  were  taken  12  hours 
earlier  with  solar  noon  at  the  Greenwich  meridian,  giving  an 
hour  angle  for  line  4  of  35° . 
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EX 


The  composition  of  solar  insolation  was  partitioned  in 
this  study,  in  accordance  with  Joseph  (1971),  by  dividing 
the  insolation,  F(2)  at  level  2,  into  two  parts,  one  of 
which  was  subject  to  water  vapor  absorption  but  not  to  Ray- 
leigh  scattering,  F(A),  and  one  which  was  only  subject  to 
Rayleigh  scattering,  F(S),  in  a  clear  moist  atmosphere. 
These  partitions  may  be  easily  computed  as  follows: 

F(A)  =  0.349  F(2)  (4-4) 

F(S)  =  0.651  F(2)  (4-5) 

Cloud-reflectivity  is  inherent  in  both  F(A)  and  F(S)  solar 
energy  portions  with  the  presence  of  two  cloud  decks. 

To  further  define  the  solar  insolation  partitions,  the 
F(A)  energy  consists  of  those  wavelengths  X  >  0.9  ym  where 
absorption  by  water  vapor  and  carbon  dioxide  bands  is  the 
dominant  attenuation  process  for  clear  air.   The  F(S)  energy 
is  composed  of  all  wavelengths  1  <  0.9  yra  where  water  vapor 
absorption  is  considered  negligible. 

B.   DISPOSITION  OF  F(S)  INSOLATION 

As  described  by  Jenks  (1974),  there  are  two  reflective- 
type  parameters  which  must  be  considered  when  studying  the 
energy  contained  in  the  shorter  wavelengths,  F(S).   The 
first  parameter  is  the  Rayleigh  clear-sky  albedo,  (after 
Coulson,  1959),  which  was  empirically  formulated  using  the 
least  squares  best-fit  technique  by  Joseph  (1971)  to  be 
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a(R)  =  .085  +  .25074  [ logx 0  (|~  Sec  Z)] 

0  (4-6) 

P0  =  1013.25  mb 

Since  mean  sea  level  is  very  nearly  1000  mb ,  this  study  as- 
sumed that  tt/Po-1.   The  value  of  Sec.  Z  is  the  inverse  of 
Eq.  (4-2). 

The  second  reflective  parameter  considered  is  the  surface 
albedo,  which  is  a  function  of  the  cosine  at  the  zenith 
angle  over  ocean  surfaces.   The  surface  albedo,  as  utilized 
by  this  study,  follows  Gates  et  al.  (1971). 

a(G)  =  max  {.06,  .06  +  .54  (.7  -  Cos  Z)}  (4-7) 

As  described  in  Sec.  Ill,  four  distinct  cloud  cases  are  pos- 
sible with  a  two-layer  cloud  model.  The  disposition  of  F(S) 
under  each  of  these  cases  will  be  discussed  in  the  remainder 
of  this  subsection. 

1 .   Clear  Sky  Case 

With  clear  sky  the  F(S)  insolation  is  modified  only 
by  Rayleigh  sky  reflectivity,  a(R),  and  surface  reflectivity, 
a(G),  as  given  by  Eqs.  (4-6)  and  (4-7).   Jenks  (1974)points 
out  that  with  multiple  reflections  possible  between  the 
earth  and  atmosphere,  the  clear  sky  insolation  absorbed  at 
the  earth's  surface  (k=10)  after  scattering  is  given  by 

IS10(0,0)  =  F(S)[l-a(R)][l-a(G)]/[l-a(R)a(G)]       (4-8) 

where  the  notation  (0,0)  denotes  CL(1)  =  CL(2)  =  0.0  as  dis- 
cussed in  section  III. 
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2.   Cloudy  Sky  Cases 

With  the  presence  of  cloud  layers,  an  additional  in- 
solation modifier  must  be  considered,  cloud  reflectivity. 
In  this  study  the  values  suggested  by  CD.  Rodgers  (1967) 
were  adopted,  namely,  R(l)  =  0.54  for  the  higher  clouds  (be- 
tween k=4  and  k=6)  and  R(2)  =  0.66  for  the  lower  clouds  (be- 
tween k=8  and  k=9). 

Consider  first  the  case  (1,0)  where  the  fraction  of 
the  grid  area  has  overcast  clouds  only  at  the  high  level, 
that  is  where  CL(1)  =  1.0,  CL(2)  =  0.0,  and  R(2)  =  0.0  so 
that  the  earth-absorbed  F(S)  insolation  is 

IS10(1,0)  =  F(S)(l-R(l))(l-a(G))/(l-R(l)a(G))       (4-9) 

The  F(S)  insolation  absorbed  by  the  surface  above 
which  there  is  present  only  a  lower  cloud  overcast,  [denoted 
by  (0,1)]  where  CL(1)  =  0.0,  CL(2)  =  1.0,  and  R(l)  =  0.0  is 
given  by 

IS10(0,1)  =  F(S)(l-R(2))(l-a(G))/(l-R(2)ct(G))      (4-10) 

If  overcast  clouds  are  present  in  both  layers,  then 
the  conditions  CL(1)  =  CL(2)  =1.0  [or  in  abbreviated  symbol- 
ism, case  (1,1)]  applied,  and  reflectance  also  exists  between 
the  cloud  layers  as  well  as  between  clouds  and  space.   Thus 
the  absorbed  F(S)  insolation  at  the  earth's  surface  as  for- 
mulated after  Arakawa  (1972)  is  given  by 

IS10(1,1)  =  F(S)[(l-R(l))(l-R(2))(l-a(G))]/[l-((R(l)R(2) 

+R(2)a(G)+R(l)a(G))+2('R(l)R(2)a(G))]     (4-11) 
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a  formula  which  resembles  very  closely  a  combination  of  Eqs. 
(4-9)  and  (4-10). 

Thus  far  in  this  subsection  the  sky  was  considered 
either  overcast  or  clear  at  each  level  without  regard  for 
the  actual  cloud  conditions  CL(1)  and  CL(2). 
3.   Composite  F(S)  Insolation 

If  the  areal  coverage  weights,  as  defined  by  Eq. 
(3-4),  are  applied  to  the  appropriate  right  sides  of  Eqs. 
(4-8),  (4-9),  (4-10),  and  (4-11),  then  the  composite  surface- 
absorbed  insolation,  IS10,  results 

IS10  =  IS10(0,0)W(0,0)+IS10(1,0)W(1,0) 

(4-12) 
+  IS10(0,1)W(0,1)+IS10(1,1)W(1,1) 

It  is  apparent  that  F(S)  insolation  reflected  to  space  may 
be  found  by  subtracting  IS10  from  F(S). 

C.   DISPOSITION  OF  F(A)  INSOLATION 

This  subsection  deals  with  that  portion  of  the  solar  in- 
solation which  is  absorbed  by  atmospheric  water  vapor  and 
carbon  dioxide,  namely  the  longer  wavelengths.   The  discus- 
sion will  again  be  divided  into  the  four  different  cloud- 
layer  cases. 

1.   The  Clear  Sky  Case 

Formulation  enabling  the  determination  of  absorption 
values  for  particular  layers  of  interest  was  provided  by  use 
of  the  Manabe-Moller  absorptivity  function  which  may  be  gen- 
eralized as  follows 

a(2,2k)  =  .271[U(2,2k)  Sec  Z]'303  (4-13) 
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Here  absorptivity  a  is  applied  over  the  pressure-scaled  water 
vapor  mass  between  levels  2  and  2k  (Fig.  1)  along  the  slant 
path  of  the  zenith  angle  Z.   Recalling  that  the  water  vapor 
mass  above  level  2  was  assumed  negligible,  the  absorbed  solar 
insolation,  for  example,  in  the  layer  (2,6)  is  written  as 

A(2,6)  =  0.271  F(A)[U(2,6)  Sec  Z]-303  (4-14) 

and  the  absorbed  solar  insolation  in  the  layer  (6,10)  is 
given  by 

A(6,10)  =  A(2,10)  -  A(2,6)  (4-15) 

For  a  schematic  depiction  of  the  F(A)  insolation  disposition 
in  the  clear  sky  case  [see  Fig.  3  in  Jenks  (1974)] . 

To  obtain  that  portion  of  F(A)  insolation  which  is 
absorbed  at  the  earth's  surface,  the  atmospheric  transmissiv- 
ity  (l-a(2,10))  and  the  surface  transmissivity  after  reflec- 
tance are  combined  (by  the  product-law  of  transmissivities) 
to  give  IA10  as  follows: 

IA10(0,0)  =  F(A){1-.271[U(2,10)  Sec  Z] • 3 ° 3 } ( l-a(G) ) 

(4-16) 

Here  a(G)  is  defined  by  Eq.  (4-7) 

In  the  remainder  of  this  subsection,  which  discusses 
the  cloudy  layer  cases,  representative  cloud  reflectivities 
and  cloud  absorptivit ies  were  adopted,  after  CD.  Rodgers 
(1967),  for  the  two  possible  cloud  layers.   The  cloud  re- 
flectivities were  RA(1)  =  .46  and  RA(2)  =  .50  while  the 
cloud  absorptivities  were  A(l)  =  .20  and  A(2)  =  .30.   Note 
that  the  reflectivities  for  the  F(A)  wavelengths  differ  from 
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those  adopted  for  the  F(S)  wavelengths.   The  procedure  of 
considering  the  cloud  conditions  to  be  overcast  whenever 
they  appear  and  then  applying  the  appropriate  weighting  fac- 
tors in  the  composite  summation  will  again  be  followed  as  in 
Sec.  IV. B. 

To  simplify  the  discussion  for  the  cloud-covered 
cases,  the  (1,1)  case  will  be  presented  first  as  it  contains 
representative  type  equations  for  the  remaining  two  cases, 
(1,0)  and  (0,1). 

2.   Overcast  in  Both  High  and  Low  Cloud  Layers 

The  following  set  of  formulas  detail  the  disposition 
of  incoming  solar  insolation  (F(A))  from  level  k=2  to  the 
earth's  surface  to  permit  determination  of  the  amount  of  in- 
solation absorbed  by  the  atmospheric  layers  and  by  the  earth's 
surface.   The  dashed  separation  lines  are  introduced  to  sub- 
divide the  absorption  and  reflection  physics  of  the  model 
into  subsections  which  permit  the  analysis  to  proceed  more 
or  less  within  successive  200  mb  layers  (see  also  Fig.  3). 

F44-  =  F(A)(1-.271[U(2,4)  Sec  Z]*303) 
F4  +  =  F4!(RA(1)) 

A244-  =  F(A).271[U(2,4)  Sec  Z]-303 
A24+  =  F4+.271[U(2,4)  Sec  Z]*303 

F6+  =  F4+(1-RA(1)-A(1)) 

A46  =  F4+(A(1))  (4-17b) 

A26  =  A244-  +  A24  +  +  A46 


(4-17a) 
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Figure  3.   Schematic  representation  of  F(A)  insolation 
disposition  in  the  case  of  two  overcast  layers. 
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TD68    =    1-. 271[U(6, 8)5/3]* 303 
TD910   =    1-. 271[U(9, 10)5/3]  '  303 


(4-18) 


F84-    =    F6|(TD68) 

F8+    =    F8+(RA(2)) 

F6+    =   F8+ (TD68) 

F64-4-    =    F6+(RA(1)) 

F8ii    =   F64-+(TD68)  (4-19a) 

A68I    =   F64-    -    F84- 

A68+    =    F8+    -    F6+ 

A684-4-    =    F6II    -    F84-  + 

A68    =   A68I    +   A68i    +   A68I* 


F9+    =   F8i(l-RA(2)    -    A(2))    +   F8+4- (1-A(2) ) 

(4-19b) 

A89   =    F8+(A(2))    +   F8++(A(2)) 


F1CH    =    F9I(TD910) 

FlOt    =   Fl(H(a(G)) 

F9t    =   F10+(TD910) 

F94-4-    =    F9+(RA(2)) 

F104-+    =   F94-KTD910)  (4-19c) 

A91CH    =   F9+    -    F104- 

A910t  =  FlOt  -  F9+ 

A910I4-  =  F94-4-  -  F104-I 

A910  =  A910+  +  A910+  +  A91CH4- 

A810  =  A89  +  A910 
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In  the  preceding  equations  the  symbol  "A"  followed 
by  a  set  of  digits  refers  to  the  F(A)  insolation  absorbed  in 
the  indicated  layer,  for  example,  A26  indicates  the  absorp- 
tion in  the  layer  k=2  to  k=6.   The  symbol  "F"  followed  by  a 
k  value  indicates  the  amount  of  insolation  passing  through 
the  particular  k  level.   The  vertical  arrow  modifiers  indi- 
cate the  direction  of  insolation  passage,  i.e.,  4-  downward 
insolation,  f  upward  reflected  insolation,  and  4-1  downward 
reflected  insolation.   Terms  involving  the  symbols  "TD"  as 
expressed  by  the  functions  of  Eq.  (4-18),  indicate  the 
Manabe-Mbller  transmissivities  for  insolation  beneath  a 
cloud  layer.   In  the  latter  situation,  the  term  Sec  Z  in 
Eq.  (4-13)  is  effectively  replaced  by  the  mean  slant-path 


Sec  Z  -  5/3  (Katayama,  1966). 

Note  that  the  effect  of  multiple  reflections  between 
clouds  or  between  the  earth's  surface  and  the  lower  cloud 
has  been  incorporated  to  the  extent  of  two  reflections,  with 
the  lowermost  reflecting  surface  absorbing  the  remaining  im- 
pinging insolation.   Computations  indicated  that  the  insola- 
tion remaining  after  two  reflections  was  too  small  to  warrant 
further  reflections.   Also  insolation  reflected  upward  from 
a  lower  interface  (cloud  or  ground)  to  the  base  of  an  upper 
cloud  deck  has  not  been  subjected  to  cloud  absorption.   This 
tends  to  reduce  very  slightly  the  secondary  cloud-absorption. 

Note  finally  that  between  any  two  reflecting  surfaces 
there  are  three  contributions  to  absorption  of  solar  insola- 
tion within  the  layer.   These  contributions  are  denoted  by 
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the  arrow-identifiers  which  indicate  the  segment  of  path 
being  traversed. 

From  Eqs.  (4-19c)  the  impinging  F(A)  insolation  at 
the  earth's  surface  can  be  expressed 

TRANA(1,1)  =  F10+  +  F10+  +  (4-20) 

The  F(A)  insolation  which  is  actually  absorbed  by  the  earth's 
surface  (see  Fig.  3)  may  be  written  as 

IA10(1,1)  =  Fl(H(l-a(G))  +  F10+  +  (4-21) 

3 .   Disposition  of  F(A)  Insolation  with  an  Upper  Overcast 
With  a  cloud  layer  present  only  in  the  upper  layer 
the  equations  depicting  the  disposition  of  incoming  insola- 
tion becomes  a  simplified  subset  of  the  previous  case: 

A24+  =  F(A).271[U(2,4)  Sec  Z]'303 

A24+  =  F(A){1-.271[U(2,4)  Sec  Z]'303}RA(1)        (4-22a) 

x{ ,271[U(2,4)  Sec  Z] ' 303} 
A24  =  A244-  +  A24t 

F6i  =  F(A){1-.271[U(2,4)  Sec  Z] • 3 ° 3 } [ 1-RA( 1)-A( 1 ) ] 

A46  =  F(A){1-.271[U(2,4)  Sec  Z]*303}A(1)  (4-22b) 

A26  =  A24  +  A46 
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F1CH    =    F6KTD610) 
FlOt  =   Fl(H(a(G)) 

F6  +    =    F10t(TD610) 

F6I+   =    F6+(RA(1)) 

F1CH  +  =  F6H(TD610)  (4-22c) 

A61CH  =  F64-  -  F1CH 

A610  +  =  F1CH  -  F6+ 

A61CH+  =  F64-I  -  F104-4- 

A610  =  A610  +  +  A61CH  +  A61(m 

The  variables  used  above  are  defined  in  a  similar  manner  to 
those  in  the  (1,1)  case.   The  impinging  insolation  at  the 
earth's  surface  may  be  formulated  as  follows 

TRANA(1,0)  =  F104-  +  F10I4-  (4-23) 

while  the  F(A)  insolation  absorbed  by  the  earth  in  this  case 
(Fig.  4)  is  given  by 

IA10(1,0)  =  Fl(H(l-a(G))  +  F10I+  (4-24) 

4.   Disposition  of  F(A)  Insolation  with  a  Low  Overcast 
With  an  overcast  lower  cloud  layer  the  disposition 
formulas  are  as  follows: 
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Figure  4.   Schematic  representation  at  F(A)  insolation  dis- 
position with  an  upper  overcast  layer  only. 
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A26I    =   F(A).271[U(2,6)    Sec   Z]  •  3  °  3 

A68I    =   F(A).271{[  U(2,8)    Sec   Z]  "  3  °  3  -[  U(2  ,  6)    Sec   Z]  •  3  °  3 } 

F8+    =   F(A){  1-.271[  U(2,8)    Sec   Z]  •  3  ° 3  } 

F8t    =    F8|[RA(2)1 

F6+    =   F8t{l-.271[U(6,8)    Sec   Z]  •  3 ° 3 }  (4-25a) 

F2+    =   F8+{1-.271[U(2,8)    Sec   Z]  •  3  ° 3 } 

A68+    =   F8+    -    F6t 

A26t    =   F6t    -    F2+ 

A26   =   A264-    +   A26+ 


F9I    =    F8I[1-RA(2)    -    A(2)] 
A89    =    F8+[A(2)] 


(4-25b) 


F10I    =    F9+(TD910) 

FlOt    =   F10l(a(G)) 

F9t    =    F10t(TD910) 

F94-+    =    F9+(RA(2)) 

F1CH  +    =   F9I+(TD910)  (4-25c) 

A910+    =    F94-    -    F10+ 

A91CH  =  F10+  -  F9+ 

A9104-i  =  F9+  +  -  F1CHI 

A910  =  A910+  +  A910+  +  A9104-  + 

A610  =  A68I  +  A68  +  +  A89  +  A910 

Again  the  variables  are  defined  as  stated  before  in  the 
(1,1)  case.   Likewise  the  incident  flux  is  defined  at  the 
earth's  surface  as 

TRANA(0,1)  =  F10+  +  FKH4-  (4-26) 
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while  that  portion  which  is  absorbed  by  the  earth,  (Fig.  5) 
is 

IA10(0,1)  =  Fl(M-(l-a(G))  +  F10II  (4-27) 

5.  Weighted  F(A)  Layer-Absorptions  and  Surface-Absorp- 
tion Insolation 

As  described  in  Section  IV. B. 3  on  the  composite  F(S) 
insolation,  the  values  of  the  absorbed  F(A)  insolation  in 
the  layers  and  of  the  earth's  surface  which  were  calculated 
for  the  four  cloud  cases  (0,0),  (1,0),  (0,1),  and  (1,1)  must 
be  appropriately  weighted  by  the  values  from  Eq.  (3-4)  to 
obtain  the  composite  F(A)  layer-  and  surface-absorbed  inso- 
lation . 

A26  =  A26(Q  Q)W(0,0)+A26(1  Q)W( 1 , 0)+A26( Q ^ 1 }W( 0 , 1 ) 

(4-28) 
+A26(11)W(1,1) 

A610  =  A610(Q  Q)W(0,0)+A610(1  Q)W( 1 , 0)+A610( Q    1)W(0,1) 

(4-29) 
+A610(1  1)W(1,1) 

IA10  =  IA10(0,0)W(0,0)+IA10(1,0)W(1,0)+IA10(0,1)W(0,1) 

(4-30) 
+IA10(1,1)W(1,1) 

6.  Absorptivity  (ABA)  by  Layers 

It  was  convenient  to  work  with  the  absorptivity  as 
well  as  the  actual  insolation  energies  absorbed  by  each  lay- 
er and  at  the  surface.   To  obtain  the  absorptivity,  the 
total  incoming  insolation  at  k=0  must  be  known.   This  value 
called  FADJ  is  formulated  as  follows 
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Figure  5.   Schematic  representation  of  F(A)  insolation 
disposition  with  a  lower  overcast  layer  only. 
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FADJ  =  2.00  (r/rm)-2  Cos  Z  (4-31) 

By  taking  the  ratio  of  the  insolation  absorbed  in  the  tropo- 
sphere to  that  incident  at  the  top  of  the  atmosphere,  the 
full  tropospheric  absorptivity  is  formed 

ABA  .  M6_±_A610  (4.32) 

Similar  formulas  may  be  expressed  for  the  partial  absorptiv- 
ities  for  the  layers  (2,6)  and  (6,10). 

D.   REFLECTION  OF  F(A)  AND  F(S)  ENERGY 

To  determine  the  albedo  and  fractional  absorptivity  of 
the  earth's  surface,  both  the  F(A)  and  F(S)  fractions  of  so- 
lar insolation  must  be  considered. 

1.   Albedo 

The  amount  of  F(A)  insolation  which  is  reflected 
back  to  space  is  determined  by  the  following 

REFA  =  F(A)  -  A26  -  A610  -  IA10  •  (4-33) 

here  A26,  A610,  and  IA10  are  defined  by  Eqs.  (4-28),  (4-29), 
and  (4-30),  respectively. 

Because  the  F(S)  wavelengths  do  not  undergo  atmos- 
pheric absorption,  the  reflected  F(S)  insolation  was  easier 
to  calculate 

REFS  =  F(S)  -  IS10  (4-34) 

IS10  was  defined  by  Eq .  (4-12).  By  combining  Eqs.  (4-33) 
and  (4-34)  the  total  reflected  insolation  to  space  is  de- 
fined by 
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REF  -  REFA  +  REFS  (4-35) 

The  resulting  REF  values  were  calculated  at  each  gridpoint 
and  are  meridionally  displayed  by  latitudinal  interpolation 
in  Figs.  7(a) , . . . , 10(b) . 

By  taking  the  ratio  of  the  total  reflected  energy  to 
the  total  incoming  energy  the  earth-atmospheric  system  albedo 
is  defined 

ALB  =  REF/FADJ  (4-36) 

where  FADJ  is  defined  by  Eq .  (4-31). 

2.   Composite  Absorptivity  (ABG)  by  the  Earth's  Surface 
The  fractional  absorptivity  of  the  earth's  surface 
was  determined  by  the  ratio  of  the  total  insolation  absorbed 
by  the  earth  to  the  total  incoming  insolation.   The  total 
insolation  absorbed  by  the  surface  is  determined  by  the  prop- 
erly weighted  sums  of  IA10  and  IS10,  defined  by  Eqs.  (4-30) 
and  (4-12),  respectively. 

ABG  =  (IA10  +  IS10)/FADJ  (4-37) 

The  result  of  Eq .  (4-37)  gives  the  weighted  absorptivity  of 
earth.   It  is  also  of  interest  to  determine  the  total  inso- 
lation incident  at  the  earth's  surface  just  before  absorp- 
tion by  earth.   The  result  is  given  by  the  properly  weighted 
sum  of  terms  of  the  form 

TRAN  =  TRANA  +  [  IS10/(  l-ct(G) )  ]  (4-38) 
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It  should  be  noted  that  each  of  the  four  terms  of  form  IS10 
(Eq.  4-12)  had  a  factor  (l-a(G))  in  the  numerator  so  that  the 
transmitted  F(S)  insolation  incident  at  earth  before  absorp- 
tion need  only  be  divided  by  (l-a(G)).   The  parameter  TRAN 
is  a  simulation  of  the  insolational  energy  incident  at  a 
pyrheliometer .   Note  finally  that  the  ratio 

ATRAN  E  TRAN/FADJ  (4-39) 

represents  the  transmissivity  of  the  troposphere. 
3.   Computational  Check 

Recalling  the  assumption  of  a  four  percent  attenua- 
tion of  the  incoming  solar  insolation  in  the  stratosphere, 
the  summation  of  the  fractional  absorptivities  ALB,  ABA,  and 
ABG  should  equal  0.96.   This  value  provides  a  good  computa- 
tional check  for  this  model  at  each  gridpoint  and  was  in 
fact  satisfied  over  the  entire  grid. 

E.   STATISTICAL  TESTS 

To  further  verify  the  foregoing  calculations,  statisti- 
cal computations  were  performed  on  several  key  parameters 
utilizing  linear  regression  techniques  from  the  BMD  set  of 
statistical  computer  programs  (Dixon,  1973). 

1.   Clear  Sky  Case  Regressions 

Three  linear  regressions  were  computed  for  the  clear 
sky  case,  the  first  of  which  was  the  clear  sky  albedo, 
ALB(0,0,Z).   The  best-fit  equation,  using  the  secant  of  the 
zenith  angle,  Sec  Z,  as  the  predictor,  was  computed  to  be 

ALB(0,0,Z)  =  0.00G50  +  0.09052  Sec  Z  (4-40) 
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with  a  regression  coefficient  R  =  .964.   This  highly  signif- 
icant correlation  coefficient  confirms  the  model  results 
which  indicate  that  Sec  Z  is  a  good  predictor  of  clear-sky 
albedo.   In  this  connection,  it  should  be  recalled  that  a(G) 
and  a(R)  were  modeled  to  be  functionally  related  to  Sec  Z. 
The  mean  global  clear-sky  albedo  was 


ALB(0,0)  =  0.164 
which  corresponds  to  the  mean  secant  of  the  zenith  angle, 


Sec  Z  =  1.743,  for  this  study.   The  statistics  show  that  the 
largest  contributions  to  error  in  the  fit  by  Eq.  (4-40)  oc- 
curred at  largest  zenith  angles,  which  occur  at  high  lati- 
tudes. 

The  two  other  clear-sky  linear  regressions  performed 
were  ABA(0,0,M)  and  ATRAN( 0 , 0 , M) ,  where  the  predictor,  M, 
was  the  water-vapor  mass  path  length  given  by 

M  =  (U  Sec  Z)5[log10(U  Sec  Z)*]  (4-41) 

This  parameterization  of  the  water-vapor  mass  is  formulated 
after  Hanson  (1967,  1971),  who  used  a  similar  absorber-mass 
predictor  in  his  empirical  equations  for  ABA  with  both  clear 
and  partly  cloudy  skies.   The  resulting  clear-sky  regression 
equations  in  this  study  for  ABA  and  ATRAN  were 

ABA(0,0,M)  =  0.10320  +  0.05535M  (4-42) 

with  R  =  .9925  and 

ATRAN(0,0,M)  =  0.77494  -  0.04823M  (4-43) 
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with  R  =  .7142.   The  sample  means  of  the  statistics  which 
entered  into  (4-42)  and  (4-43)  are 


ABA(0,0,M)  =  0.146 


ATRAN(0,0,M)  =  0.738 


M  =  0.776  (gm  cm"2 )* 


based  on  the  original  93  cases. 

2.   Cloudy-Sky  Case  Regressions 

In  addition  to  the  clear-sky  fractional  amounts  of 
ALB,  ABA,  and  ATRAN,  linear  regressions  were  formulated  for 
the  composite  cloudy-sky  case,  i.e.  ALB[ CL( 1) ,CL(2) ] , 
ABA[CL(1),CL(2)] ,  and  ATRAN[ CL( 1 ) ,CL( 2 ) ] ,  respectively.   The 
regression  relationships  between  the  clear-  and  cloudy-sky 
cases  were  determined  for  the  same  zenith  angle  at  gridpoint 
(I, J)  using  the  total  opaque  cloud  cover  (CL)  as  the  predic- 
tor.  The  total  opaque  cloud  cover  approximates  the  effective 
cloud  cover  existing  in  a  two-layer  cloud  model  and  is  ex- 
pressed as 

CL  =  CL(1)  +  CL(2)  -  CL(1)CL(2)  (4-44) 

The  relationship  between  the  cloudy-sky  albedo  and 
the  corresponding  clear-sky  calculation  was  obtained  by 
least  squares  with  the  following  results 

ALB[CL(1) ,CL(2)]  =  ALB( 0 , 0) [ 1+2 . 09929CL+0. 84608CL2 1    (4-45) 

and  R  =  .9292.   The  means  of  the  parameters  used  in  Eq.  (4-45) 
are 
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ALB[CL(1),CL(2)]  =  0.406 


ALB(0,0)  =  0.164 
CL  =  0.662 

The  dependence  of  Eq.  (4-45)  upon  CL  is  indicated  by 
the  high  correlation  coefficient  (R=.9292).   The  global  mean 
albedo  for  this  study  of  0.406  is  somewhat  higher  than  that 
reported  in  recent  literature  (ALB=.300)  by  Raschke  et  al. 
(1973).   However,  Raschke' s  figure  is  an  annual  mean  which 
includes  the  somewhat  smaller  cloud  amounts  over  land. 

The  tropospheric  absorptivity  of  air  in  the  composite 
cloudy  case  was  developed  analogously  to  (4-45)  by  least 
squares  and  gave  the  form 

ABA[CL(1) ,CL(2)]  =  ABA( 0 , 0) [ 1+0 . 22305CL+0 . 17261CL2 ]    (4-46) 

with  R  =  .927  and  mean  statistics  for  the  above  parameters 


ABA[CL(1),CL(2)]  =  0.180 


ABA(0,0)  =  0.146 

CL  =  0.662 

The  direct  relationship  between  solar  absorptivity  and  cloud 
cover  is  in  general  agreement  with  the  results  of  Plante 
(1972)  who  found  increased  absorptivity  with  increasing  CL. 
This  relationship  could  be  the  reason  for  the  recent  reduced 
estimates  of  global  albedo  [Raschke  et  al.  (1973),  Vonder 
Haar  and  Suomi  (1971)].   The  significant  value  of  R  indi- 
cates the  degree  of  sensitivity  between  cloud  amounts  and 
atmospheric  absorptivity. 
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The  final  statistical  regression  equation  formulated 
deals  with  the  cloud-sky  transmissivity  of  the  troposphere 
relative  to  the  clear-sky  calculation.   Again  the  total 
opaque  cloud  cover  has  been  used  as  the  independent  variable 

ATRAN[CL(1),CL(2)]  =  ATRAN( 0 , 0) [ 1-0 . 51477CL-0 . 14646CL2 ] 

(4-47) 

Here  the  correlation  coefficient  R  =  .997  and  the  mean  sta- 
tistics over  the  93  sample  soundings  are 


ATRAN[CL(1),CL(2)]  =  0.543 


ATRAN(0,0)  =  0.738 

CL  =  0.662 

Note  the  high  degree  of  significance,  (R  =  .997)  of  Eq. 
(4-47),  which  formulation  agrees  with  empirical  results  of 
Quinn  (1971)  and  numerous  other  workers  cited  by  Quinn  in 
his  1971  paper. 

In  general,  the  solar-model  statistical  relation- 
ships developed  in  this  study  show  consistency  with  empiri- 
cally reported  observations,  which  have  already  been  noted, 
when  consideration  is  taken  for  seasonal  effects  and  the 
latitudinal  limits  of  this  study. 
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V.   SENSIBLE  AND  LATENT  HEAT  TRANSPORT 
AT  THE  SEA- AIR  INTERFACE 

A.  GENERAL  CONSIDERATIONS 

For  a  complete  discussion  of  the  heat  balances  of  the 
earth  and  atmosphere  for  16  October  1973,  the  vertical 
transports  of  sensible  and  of  latent  heats  across  the  sea- 
air  interface  must  be  considered.   For  this  purpose,  the 
turbulent  flux  models  for  the  planetary  boundary  layer,  as 
used  in  the  FNWC  primitive  equation  model   as  described  by 
Kesel  and  Winninghoff  (1972)  and  Kaitala  (1974)  ,  was  adapted 
here  virtually  without  change. 

B.  EVAPORATION 

The  determination  of  the  surface  latent  heat  flux,  E, 
over  the  ocean  is  based  upon  its  large  scale  parameteriza- 
tion, principally  after  Langlois  and  Kwok  (1969) 

E  =  Lp10 K*  _9     [^  "'  ]  (5-1) 

Zj9-Zi  i  o 

Here  L  is  the  latent  heat  of  vaporization  given  by 

L  =  596.73  -  .601  Tio    (cal/gm) 
The  other  parameters  in  (5-1)  are  formulated  as  follows: 

Pio  -  Pio/RdT10 

1000. ?h 
°9    l9(  900) 


es  =  t10 
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Zo-Z 


i  o 


Rd 

9— 10           1000 
T           ln      900 

g 

691 

e 

Q7    —    -nfQr\Pt\ 

qx-q9  =  .62197  3-^  -q(900) 

Note  that  at  level  k=10,  the  assumption  P1  0  a  1000  mb  is  con- 
sidered valid  and  e   is  the  vapor  pressure  at  the  top  of  the 
constant  flux  layer.   Also  q(900)  is  the  mixing  ratio  at  k=9 
as  defined  by  Eq.  (2-4).   In  Eq.  (5-1)  K*  and  a*  are  turbulent 
transfer  coefficients  having  values  2  x  105  cm2  sec-1  and 
5  x  104  cm/°K,  respectively,  according  to  Kaitala  (1974). 
However,  values  of  the  eddy  coefficient,  K* ,  of  2  x  105  and 
1  x  105  were  both  tested  in  this  study  and  it  was  decided 
primarily  to  use  only  the  latter  value.   Furthermore,  it  is 
felt  by  this  author  that  upon  further  testing,  with  a  broader 
data  base,  it  may  be  feasible  to  tune  K*  between  the  two 
values  tested  so  as  to  make  it  dependent  either  upon  latitude 
or  season,  or  other  pertinent  variables. 

C.   SENSIBLE  HEAT  TRANSPORT 

Sensible  heating  is  computed  as  a  function  of  the  tem- 
perature difference  between  the  surface,  T10,  and  the  air 
near  the  surface,  T  ,  which  is  at  the  top  of  the  constant 
flux  layer.   Since  wind  values  were  not  carried  in  this 
study,  the  technique  used  by  Mintz  and  Arakawa,  as  described 
by  Kaitala  (1974)  and  by  Kesel  and  Winninghoff  (1972),  was 
first  attempted  for  sensible  heating 
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k*         e_-e 


Hr  =  P-V ^~mc  -  zf^V         <5-2> 

l+a*(^5— ^  ) 

Here  C   =  .239  cal/(gm°K),  and,  since  the  constant  flux  lay- 

er  was  assumed  to  be  quite  thin,  0  =T  .  y   is  a  critical 

x   x   '  c 

lapse  rate  (3  x  10~5  °K/cm),  and  the  remainder  of  the  terms 

are  defined  as  in  the  discussion  following  Eq .  (5-1).   The 

micrometeorological  equivalent  of  Eq .  (5-2)  is  given  by 
Kaitala  (1974)  as 

Hr  '■  PiiCpWVV  (5"3) 

It  can  be  seen  from  Eqs.  (5-2),  (5-3)  that  the  temperature 

lapse,  T  -T  ,  for  .the  thin  layer  between  the  surface  and  the 
g   x 

top  of  the  constant  flux  layer  becomes  very  critical  in  de- 
termining the  amount  of  sensible  heating  (cooling)  for  each 
sounding  location. 

In  Eq.  (5-2),  the  expression  within  the  brace  determined 
the  sign  of  the  sensible  heat  flux  in  the  vertical.   In  ap- 
proximately one  out  of  three  cases,  the  braced  expression 
gave  rise  to  negative  values  of  Hp ,  whereas  the  correspond- 
ing latent  heat  flux,  E,  of  Eq.  (5-1)  was  systematically 
positive.   Based  on  Fig. 15  in  Malkus  ( 1962),  which  shows  Kp>0, 
there  seems  to  be  no  valid  basis  for  having  an  erratically 
downward  distribution  of  sensible  heat  over  the  open  oceans. 
Consequently,  this  study  assumed  that  E  was  properly  as- 
sessed by  Eq.  (5-1),  and  that  a  reasonable  set  of  values  of 
Hp  could  be  obtained  by  means  of 
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Hp  =  rE  (5-4) 

Here  the  Bowen  ratio  r  was  determined  from  a  set  of  best  es- 
timates as  a  function  of  latitude. 

Using  20  values  of  r  and  their  functionally  related  lati- 
tudes, taken  from  published  results  of  Sellers  (1965), 
Budyko  (1956),  and  Jacobs  (1951),  a  linear  regression  was 
best-fitted  to  |  Sin  <J>  |  to  give 

r  =  0.10910  -  .41996  |Sin  <J>  |  +  .89645  |  Sin  <J>  |  2     (5-5) 

with  a  correlation  coefficient  of  R  =  .982.   Thus,  using 
Eqs.  (5-4)  and  (5-5),  the  sensible  heat  transport  was  ob- 
tained at  each  gridpoint,  while  always  maintaining  proper 
sign  sense  with  E  which  was  virtually  always  positive. 

D.   COMBINED  TURBULENT  HEAT  FLUX  ACROSS  THE  SEA-AIR  INTER- 
FACE 

To  get  a  better  feel  for  the  effect  that  sensible  heat- 
ing has  at  the  surface,  the  summed  terms  E  and  Hp  were  uti- 
lized as  heat-loss  quantities  in  the  total  budget  at  the 
surface,  therefore  in  terms  of  large  scale  "convective"  dis- 
position 

E  +  Hp  =  E(l+r)  (5-6) 

where  r  ranges  in  magnitude  from  approximately  .05  to  .5 
with  increasing  latitude. 
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VI .   MERIDIONAL  CROSS-SECTIONAL  DEPICTION  OF 
THE  HEATING  MODEL  COMPUTATIONS 

A.   GENERAL  DESIGN  OF  HEAT-MODEL  OUTPUT 

In  this  section,  the  computational  subroutines  already 
developed  (Sees.  Ill,  IV,  and  V)  are  combined  to  perform  a 
single  time-step  of  the  heating  program  developed  for  adapta- 
tion to  the  FNWC  prediction  model.   The  results  represent 
computations  of  the  relevant  heating  functions  at  each  grid- 
point  of  the  four  meridians  (Fig.  2). 

For  the  three  Pacific  cross-sections,  the  FNV/C  gridpoint- 
processed  analyses  for  OOOOGMT,  16  October  1973,  have  been 
used,  whereas  for  the  35° W  (Atlantic)  meridian,  the  FNWC 
analyses  from  1200GMT,  15  October  1973,  were  utilized.   For 
simplicity  the  heating  model  computations  will  be  referred 
to  as  those  pertaining  to  16  October  1973. 

Fig.  6  depicts  the  key  to  the  computations  made  at  (I, J) 
gridpoint  soundings,  shown  in  the  form  of  Table  1(b),  and  in 
effect  represents  the  flow  diagram  of  computations  proceed- 
ing from  top-to-bottom  at  each  gridpoint.   Figs.  7,  8,  9, 
and  10  show  that  the  gridpoints  considered  are  at  5°  lati- 
tude increments,  but  this  has  been  done  merely  for  simplifi- 
cation of  the  interpretations  and  for  comparisons  with 
climatological  data.   In  order  to  deduce  a  heating  parameter, 
Q(I),  at  a  non-integral  grid-value  I  corresponding  to  lati- 
tude <$>    (in  whole  multiples  of  5°),  a  Lagrangian  cubic  i  inter- 
polation scheme  was  utilized  of  the  form 
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on,    =  Qi  (I-2)(I-3)(I-4)    Qg(I-l)(I-3)(I-4) 
^XJ  (1-2)  (1-3)  (1-4)      (2-1)  (2-3)  (2-4) 

(6-1) 
Qa(I-l)(I-2)(I-4)    Q„(I-l)(I-2)(I-3) 
(3-1)  (3-2) (3-4)      (4-1) (4-2) (4-3) 

Here  I  is  determined  by  an  equation  of  the  form  (4-3c),  the 
points  1,2,3,  and  4  are  any  four  successive  integral  values 
of  I  with  I  lying  in  general  between  points  2  and  3,  although 
the  subroutine  also  works  well  if  I  is  non-centered  in  the 
I-array . 

B.   CROSS-SECTIONAL  DATA  AT  LEVEL  k=2 

The  terms  explained  here  are  those  occurring  at  the  lev- 
el k=2.   The  discussion  of  all  insolational  parameters  dis- 
cussed previously  in  section  IV  pertained  to  the  specific 
time  of  day  corresponding  to  the  hour  angle  h  at  the  instan- 
taneous times  under  consideration.   Thus  the  incident  inso- 
lation was  given  by  Eq.  (4-1),  i.e. 

F(2)  =  S  (|-)"2  Cos  Z 
m 

with  Cos  Z  given  by  Eq.  (4-2).   For  the  specific  map  times 
the  local  solar  times  gave  h=35° ,  10° ,  55° ,  and  35°  for 
cross-sections  1,2,3,  and  4  depicted  in  Figs.  7,  8,  9,  and 
10  respectively. 

For  each  specific  zenith  angle  Z  corresponding  to  lati- 
tude <j>(I,J,t)  the  absorbed  insolations  A(2,6,t)  and  A(6,10,t) 
of  Eqs.  (4-28)  and  (4-29)  have  been  determined  for  the  given 
time  t  of  16  October  1973  for  which  the  specific  hour  angles 
were  applicable.   The  instantaneous  parameter  REF(t)  of  Eq. 

(4-35)  was  determined  by  the  equation 
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REF(t)  =  F(2)  -  A26  -  A610  -  ( IA10  +  ISIO)         (6-2) 

where  the  expression  in  parentheses  is  identical  to  the  in- 
stantaneous absorbed  insolation  of  earth  given  by  the  numer- 
ator of  (4-37) . 

The  net  terrestrial  flux  F2 *  was  derived  in  Eq.  (3-11) 
and  was  not  subject  to  any  systematic  diurnal  variation. 
Thus  F2 *  may  be  taken  as  constant  for  a  24-hour  period. 

Note  that  the  first  symbol  at  k=2  in  Fig.  6  is  not  F(2), 
but  QAVE,  defined  for  each  gridpoint  by  the  analog  of  (4-1), 
namely 


r  s- 


QAVE  =  S(  — )-2  Cos  Z  (6-3) 

rm 


Here  Cos  Z  is  the  24-hour  average  of  Cos  Z,  and  is  given  by 


Cos  Z  =  [H  Sin  <J)  Sin  6  +  Cos  $  Cos  6  Sin  H]/tt      (6-4) 

where 

H  =  arc  Cos  [-Tan  <J>  Tan  6]  (6-5) 

is  the  hour  angle  from  local  noon  to  sunset  at  the  gridpoint 
Thus  the  instantaneous  values  of  incident  insolation  and  of 

reflected  insolation  have  been  converted  from  F(2)  and  REF 
to  their  24-hour  averages  denoted  by 


QAVE  =  F(2)  Cos  Z/Cos  Z  (6-6) 

and 


QREF  =  REF  Cos  Z/Cos  Z  (6-7) 

respectively . 
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The  reason  for  the  conversion  in  (6-6)    and  (6-7)  is  to 
convert  the  net  incoming  insolation  at  k=2  to  values  that 
are  comparable  in  time-scale  to  those  of  F2 * .   This  device 
has  the  advantage  of  showing  whether  the  resultant  averaged 
heat-transfer  quantities  are  compatible  with  the  expected 
approximation  to  a  heat  budget  for  the  given  date.   However, 
if  it  were  desired  to  exhibit  the  expected  diurnal  variation 
of  heat-transfers  across  the  layers  (2,6),  (6,10)  and  at  the 
ground,  the  quantity  QAVE  and  all  other  solar-heating  rates 
to  be  discussed  would  be  multiplied  by  the  factor 


(Cos  Z/Cos  Z) 

assuming  that  all  the  layer  absorptivities  and  the  solar  re- 
flectivity are  constant  during  the  day  for  which  the  diurnal 
variation  of  heat  transfers  are  of  interest. 

With  the  solar  quantities  averaged  for  the  day,  one  ob- 
tains 

BALT  =  QAVE  -  (QREF  +  F2*)  (6-8) 

as  a  24-hour  average  radiational  balance  in  cal  cm-2  min-1 
for  the  tropopause. 

C.   CROSS-SECTIONAL  RADIATIVE  TRANSFERS  IN  LAYER  (2,6) 

At  this  layer,  in  Fig.  6,  the  heat-transfers  are  assumed, 
as  in  the  FNWC  heating  model,  to  be  radiative  only.   For  con- 
venient comparison  with  climatology,  the  heating  (cooling) 
rate  is  given  by 

BAL26  =  Q26  -  F26  (6-9) 
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where  Q26  is  the  daily  average  absorption  in  layer  (2,6)  and 
is  defined  relative  to  A26  by  an  equation  analogous  to  Eq. 
(6-6).   The  terrestrial  cooling  rate  F26  is  defined  by  Eq. 
(3-12a). 

Note  that  all  values  of  BAL26  were  negative  in  all  cross- 
sectional  diagrams,  Figs.  7,..., 10. 

D.   CROSS-SECTIONAL  DATA  FOR  LAYER  (6,10) 
1.   Radiative  Transfers  Only 

Neglecting  temporarily  the  turbulent  transfers  at 
the  earth-air  interface,  the  present  heating  model  computes 
the  24-hour  average  radiative  cooling  computed  by  the  differ- 
ence equation 

BAL610  =  Q610  -  F610  (6-10) 

in  a  manner  analogous  to  Eq.  (6-9).   The  parameter  BAL610 
was  negative  at  all  gridpoints  of  Figs.  7,..., 10.   In  fact 
the  ratio  R  =  F610/Q610  varied  with  latitude  as  follows 

0°     10°     20°     30°     40°     50°     60° 
R  =  1.916   2.219   2.577   2.893   3.912   4.122   4.903 

Even  when  the  instantaneous  (near-noon)  radiation  balance 
was  considered,  the  difference  remained  at  least  slightly 
negative  in  low  latitudes  corresponding  to  cooling.   In  high- 
er latitudes,  instantaneous  cooling  rates  were  never  reversed 
by  consideration  of  the  instantaneous  computation  of  the  ra- 
diation balance  in  layer  (6,10). 
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2.  Heat  Balance  at  the  Earth's  Surface 

In  Section  V.B,  the  earth-to-air  turbulent  transfer 
of  latent  heat,  E,  has  been  modeled  after  Kaitala  (1974), 
through  Eq.  (5-1).   The  sensible  heat  transfer,  IL,  ,  given  by 

Hr  =  rE 

where  r  is  the  Bowen  ratio  determined  empirically  by  Eq. 
(5-5)  as  a  function  of  latitude. 

The  sum  of  the  turbulent  heat  transfers,  H„+E,  has 
been  subtracted  from  the  average  radiation  at  the  earth's 
surface 

R  =  QABG  -  F10*  (6-11) 

to  form  the  average  warming  rate,  BALB,  at  the  surface 

BALB  =  QABG  -  F10*  -  (Hr+E)  (6-12) 

where  QABG  is  related  to  ABG  of  Eq.  (4-37)  by  an  analog  to 
Eq.  (6-6)  and  F10*  is  defined  by  Eq.  (3-10)  which  is  taken 
as  constant  for  a  24-hour  period  as  was  the  case  with  F2*. 

The  parameter  BALB  has  been  shown. in  each  cross-sec- 
tion as  a  function  first  of  gridpoint  on  each  meridian,  and 
subsequently  interpolated  to  5°  latitude  increments  as  indi- 
cated in  connection  with  Eq.  (6-1). 

3.  Sensible  Heat  Modification  of  the  Atmospheric  Columns 
In  Figs.  7,..., 10  the  sensible  heat  transfer,  Hp , 

which  crosses  the  earth-air  interface  has  been  added  to  the 
air  only  in  the  layer  (8,10),  (cf.  Kaitala,  1974).  In  this 
connection,  the  net  radiative  loss  in  the  layer  (6,10)  has 
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been  subdivided  equally  into  two  halves,  BAL68  and  BAL810 
(excluding  H  ),  equal  to  half  the  radiative  loss  BAL610. 

The  sensible  heat  transfer  H„ ,  which  has  been  allowed 
to  diffuse  only  into  layer  (8,10),  is  shown  as  item  (n)  of 
layer  (8,10)  in  Fig.  6.   In  the  numerical  computations,  EL 
is  the  smallest  of  all  the  heat-transfers  considered  in  any 
atmospheric  column  presented  in  Figs.  7,..., 10. 

E.   MERIDIONAL  CROSS-SECTIONS  OF  THE  VERTICAL  HEAT  BUDGET 
Note  that  the  single-time  set  of  heating  computations 
for  each  of  the  four  meridians,  125°W,  170°W,  145°E,  and 
35°W,  are  shown  in  Figs.  7,..., 10  respectively.   These  fig- 
ures show  the  individual  computations  of  the  warming  and 
cooling  which  are  schematically  formatted  in  Fig.  6,  together 
with  the  layers  to  which  they  contribute  in  the  heating  mod- 
el. 

Note  also  that  the  four  meridional  cross-sections  de- 
picted each  have  two  sections:  section  (a)  depicts  results 
for  a  "tropical"  section  with  gridpoints  extending  from  20°S 
(or  the  southernmost  latitude  to  which  the  meridian  in  ques- 
tion extends)  to  25°N,  while  section  (b)  of  Figs.  7,..., 10 
depicts  the  gridpoint  computations  (averaged  for  24-hours) 
over  the  latitudinal  range  30°N  to  65°N  (or  the  northernmost 
latitude  to  which  the  meridian  in  question  extends).   Recall 
the  latitudinal  boundaries  for  the  selected  meridians  were 
required  to  remain  over  the  ice-free  ocean  on  the  FNWC  grid. 

The  latitudinal  variation  of  both  the  short-  and  long- 
wave calculations,  as  observed  in  Figs.  7,...  ,10,  seem  to 
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have  uniform  gradients  across  each  meridian  in  spite  of  the 
differing  climatic  zones  represented.   Hence  it  was  reason- 
able to  formulate  a  mean  radiational  balance  model  for  16 
October  1973,  which  is  presented  in  Section  VII. 
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VII.   MEAN  BUDGETARY  BALANCE  MODEL 

By  averaging  key  parameters,  discussed  in'  Section  VI, 
across  the  four  possible  meridional  lines  at  each  5°  latitude 
increment  from  0°  to  65° N,  a  mean  heat  budget  model  was  con- 
structed which  may  be  considered  representative  of  Fall  in 
the  annual  cycle.   The  advantage  of  such  averaging  is  to  fa- 
cilitate the  comparisons  of  this  study's  results  with  other 
published  heat  budget  results  (cf.  Budyko,  1956)  even  though 
such  results  are  generally  annual  means.   Also,  by  averaging 
over  the  four  meridians  studied  here,  meridional  differences 
tend  to  be  smoothed  out,  so  that  the  zonally-averaged  heat 
budgets  of  this  study  may  be  compared  with  others  in  the  lit- 
erature. 

A.   RADIATION  MODEL  OF  THE  OCEAN- ATMOSPHERE  SYSTEM 

Fig.  11(a)  depicts  the  schematic  format  of  the  mean  ra- 
diation meridional  cross-section  (Figs.  11(b)  and  (c)). 
Turbulent  transfer  parameters  have  been  omitted  from  this 
display  and  only  Northern  Hemisphere  values  are  given. 

Note  that  the  radiational  quantities  for  the  layers 
(6,8)  and  (8,10)  of  Figs.  7,..., 10  have  been  combined  at 
identical  latitudes  in  Fig.  11  as  for  the  layer  (6,10)  and 
that  all  resultant  values  are  meridionally  averaged  where 
possible.   Since  only  the  170°W  meridian  contains  results  at 
65°N,  the  value  at  this  latitude  is  not  representative  of  the 
zonal-averaging  process  described  above. 
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To  facilitate  a  graphical  display  of  the  results  con- 
tained in  Fig.  11  and  for  ease  of  comparison  with  published 
yearly  climatological  results  (Malkus,  1962),  the  model  dis- 
played in  Fig.  11  has  been  used  as  the  basis  for  testing  the 
following  component  heat  balances:  (a)  the  earth-atmosphere 
system  balance  (using  symbols  consistent  with  Malkus  (1962)), 
(b)  the  atmospheric  radiative  balance,  and  (c)  the  earth- 
surface  balance. 

B.   RADIATION  BALANCES  OF  THE  OCEAN-TROPOSPHERE  SYSTEM 

The  first  graph,  Fig.  12,  consists  of  a  radiation  bal- 
ance composed  of  three  parameters,  R  ,  R  ,  and  R,  which  are 
defined  (after  Malkus,  1962)  as  follows:  (1)  R   is  mean  ra- 
diative  energy  gain  (loss)  rate  of  the  ocean-troposphere 


system  which  was  referred  to  as  BALT  at  k=2  in  Fig.  11(a); 

(2)  R   is  the  mean  radiative  cooling  rate  in  the  troposphere 
a 


which  is  defined  as  the  sum  of  BAL26  and  BAL610  from  Fig. 
11(a);  and  (3)  R,  the  mean  radiative  warming  (cooling)  rate 


at  the  earth's  surface  which  was  identified  as  BAL10  in  Fig. 
11(a).   The  relationship  between  the  three  radiative  param- 
eters can  be  expressed  as 

R   =  R  +  R  (7-1) 

s        a 

Table  III  lists  quantitative  comparisons  of  R  ,  R  ,  and 

S     cL 

R,  as  found  in  this  study  for  16  October  1973,  with  values 
reported  in  Malkus  (1962)  for  the  climatological  year,  and 
values  reported  by  Jenks  (1974)  for  25  April  1973.   The  lat- 
itudinal ly-  averaged  values  are  obtained  by  weighting  each 
value  by  the  cosine  of  the  latitude 
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Radiation 
Balance 


Figure  12.   Radiational  balance  at  the  trope-pause  (R  );  in 

the  atmospheric  column  between  k=2  and  k=10  (R  );  and  at  the 

a 

surface  (R)  for  16  October  1973. 
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55 
60 
65N 


R 


(a) 


R 


a 


(b) 

/  R„  /  R 


0456/. 

0358/. 

0149/. 

0060/. 

0353/. 

0620/. 

0871/. 

0944/. 

1066/ 

1240/- 

1423/- 

1325/- 

1517/- 

1112/ 


0627/ 
0608/ 
0589/ 
0513/ 
0437/ 
0361/ 
0266/ 
0171/ 
0.0  / 
0418/- 
0570/ 
0684/ 
0760/ 
/ 


031 
024 
032 
041 
032 
062 
068 
055 
044 
.003 
009 
022 
042 
106 


1231/- 
1361/- 
1463/- 
1522/- 
1713/- 
1896/- 
1853/- 
1650/- 
1770/- 
1738/- 
1759/- 
1472/- 
1528/- 
1233/ 


1235/- 
1254/- 
1292/- 
1387/- 
1368/- 
1235/- 
1178/- 
1102/- 
1159/- 
1273/- 
1330/- 
1330/- 
1311/- 
h 


(c) 

Rj  /   R2  /  R3 


.132 
.142 
.167 
.167 
.179 
.174 
.160 
.153 
.140 
.135 
.128 
.132 
.141 
.114 


1687/ 
1720/ 
1596/ 
1463/ 
1361/ 
1276/ 
0982/ 
0705/ 
0703/ 
0499/ 
0336/ 
0146/ 
0013/ 
0120/ 


/ 
2033/ 

/ 
2090/ 

/ 
1900/ 

/ 
1482/ 

/ 
1026/ 

/ 
0646/ 

/ 
/ 


154 
158 
194 
208 
206 
236 
225 
201 
175 
132 
124 
155 
184 
219 


Wt. 

Avg.  -.0552/. 0186/. 039   -. 1589/- . 1264/- . 150   . 1037/ . 1625/ . 185 

Table  III.   Comparisons  of  Radiational  Balances  (a)  at  the 

tropopause  R  ,  (b)  in  the  vertical  column  R   between  k=2  and 
s  a 

k=10,  and  (c)  at  the  earth's  surface  R  of  the  meridionally 
averaged  values  as  found  in  this  study  denoted  by  subscript 
"1".   The  corresponding  annual  values  reported  by  Malkus 
(1962)  are  denoted  by  subscript  "2",  and  those  for  25  April 
1973  found  by  Jenks  (1971)  denoted  by  subscript  "3".   All 
values  are  in  ly  min-1. 
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n 


I       Q   COS  (J) 
Wt.  Avg.  =  ~± - -  (7-2) 


.  £.,  Cos  A 
1=1      y 


It  can  be  noted  from  Table  III  that  although  the  R   val- 
ues  found  by  this  study  are  consistently  lower  than  those  re- 
ported by  Malkus  for  the  year,  the  seasonally  higher  values 
reported  by  Jenks  tend  to  offset  this  difference.   The  same 
holds  true  for  the  balance  at  the  surface  (R).   However, 
both  seasonal  studies  show  greater  radiative  cooling  in  the 

troposphere  (R  )  than  the  yearly  climatological  values  which 
a 

evidently  was  obtained  after  Budyko  (1956)  as  a  residual  be- 
tween R   and  R,  as  functions  of  latitude.   It  is  noteworthy 

s       '  J 

that  Jenks'  (1974)  R   results  (column  (c)  of  Tabic  III)  are 
greater  than  those  obtained  here  by  .094  ly  min-1  (in  the 
mean),  and  of  this  surplus  the  additional  heating  rate 
AR  =  0.081  ly  min-1  is  found  at  the  surface. 

C.   TROPOSPHERIC  HEATING 

Fig.  13  graphically  depicts  the  mean  tropospheric  heat 
budget  as  found  in  this  study.   The  parameters  presented  in- 
clude: R   which  was  discussed  in  Sec.  VII. B,  (E+Hr)  which  is 
a  i 

the  turbulent  heat  transfer  rate  due  to  latent  and  sensible 

heat  transport  realized  in  the  atmosphere,  Q   which  is  the 

^  *     '   va 

net  heat  flux  divergence  rate  in  the  troposphere  required  to 

bring  about  a  heat  balance  (which  would  not  be  expected  to 

be  achieved  for  the  time  scale  here),  and  S   the  storage 

a 

term.   These  terms  are  related  by 
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Heating  rales, 
Troposphere 


5         10         15         20        25        30        35        40       45        50        55        60        65 

Latitude  (°N) 
Figure   13.      Trope-spheric   Heat   Budget   Disposition.      The   net 

heat    loss    in    the  Troposphere    (Q      +Sr  )    is   depicted   as   the   sum 

of   the   tropospheric   radiative   cooling    (R    )    and   the   turbulent 

a 

transfer  heating  (E+Hr)  in  the  troposphere. 
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(a) 


(b)  (c) 


°LAT 
0 

R    / 

R 
a2 

(E+H  )  /(E+H  ) 

1   1         •   2 

(Q   +S  )  /(Q   +S  ) 
va   a  i    va   a  2 

-.1231/- 

.1235 

.0294/ 

-.0937/  .0874 

5 

-.1361/- 

.1254 

.0332/. 1539 

0.1029/  .0285 

10 

-.1463/- 

.1297 

.0425/ 

-.1038/  .0437 

15 

-.1522/- 

.1387 

.0636/. 1786 

-.0916/  .0095 

20 

-.1713/- 

.1368 

.0544/ 

-.1169/-. 0057 

25 

-.1896/- 

.1235 

.0744/. 1748 

-.1095/-. 0095 

30 

-.1853/- 

.1178 

.1164/ 

-.0749/  .0190 

35 

-.1650/- 

.1102 

.1701/. 1615 

-.0051/  .0190 

40 

-.1770/- 

.1159 

.1544/ 

-.0226/  .0038 

45 

-.1738/- 

.1273 

.1518/. 1254 

-.0220/-. 0095 

50 

-.1759/- 

.1330 

.1486/ 

-.0273/-. 0152 

55 

-.1472/- 

.1330 

.0880/. 0988 

-.0592/-. 0228 

60 

-.1528/- 

.1311 

.0924/ 

-.0604/-. 0266 

65N 

-.1233/ 

.0188/ 

-.1045/ 

Wt 

• 

Avg. 

-.1589/- 

.1264 

.0855/. 1533 

-.0735/  .0201 

Table  IV.   Comparisons  of  trapospheric  radiative  cooling  rate 

(R  ),  turbulent  transport  warming  rate  (E+H),  and  net  trop- 
a  r 

ospheric  heat  loss  rate  (Q   +S  )  as  found  by  this  study  for 

v  va   a  J  J 

16  October  1973  denoted  by  subscript  "1"  and  as  reported  by 
Malkus  (1962)  for  yearly  climatology  denoted  by  subscript  "2" 
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Q    +  S   =  R   +  (E+Hp)  (7-3) 

va    a    a    v    r  K         ' 

The  storage  rate  S   has  not  been  separated  specifically  from 

Q   ,  but  over  the  annual  cycle  it  is  assumed  zero.   Assuming, 

however  that  Q   =0,  if  the  weighted  average  values  from 
va 

Table  IV  are  set  into  Eq .  (7-3),  to  determine  the  mean  stor- 
age cooling  rate,  S  ,  the  latter  turns  out  to  be 

a 

S~   =    (||)   =  -.54°C  (day)"1 

for  the  time  (16  October  1973)  under  study.   The  net  lateral 
flux-divergence  term  Q   has  not  been  assessed  here,  but  pre- 
sumably its  mean  effect  would  show  up  as  cold  advection 
across  the  equator,  and  would  result  in  a  reduced  value  of 

S~~. 
a 

In  Fig.  13,  the  dashed  portion  of  the  (E+EL)  and  (Q   +S  ) 

1         va   a 

curves  are  displayed  to  show  the  effects  of  varying  the  K* 
value  in  the  calculation  of  latent  heat  transport  (Eq.  5-1), 
as  discussed  in  Section  V.   In  the  dashed  portion,  the  K* 
value  of  2  x  105  cm2  sec-1  was  tested  for  latitudes  0°,  5°, 
and  10°N,  with  a  transitional  value  of  1.5  *  105  cm2  sec-1 
used  at  15°N.   At  all  points  on  the  solid  (E+Hr)  line  K*  was 
1  x  105  cm2  sec-1.   This  smaller  K*  value  seemed  to  provide 
more  realistic  values  of  (E+Hp)  in  mid-latitudes  but  to  some- 
what unrealist ically  small  values  in  low  latitudes. 

A  possible  explanation  for  the  small  (E+H  )  values  at 
low  latitudes  is  the  difficulty  in  correctly  determining  the 
temperature  and  humidity  gradients  between  the  surface  and 
the  top  of  the  constant  flux  layer  in  the  data-sparse 
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equatorial  areas.   In  these  tropical  latitudes  some  time- 
and  space-smoothing  of  the  few  available  radiosondes  could 
have  been  introduced  into  the  original  FNWC  map  analyses, 
which  would  tend  to  smooth  out  vertical  gradients  as  against 
those  to  be  expected  with  observed  radiosondes  on  16  October 
1973.   Hence  it  was  considered  justifiable  to  revert  to 
Kaitala's  (1974)  suggested  K*  =  2  x  105  cm2  sec-1  in  low 
latitudes,  0°-10°N,  where  the  resultant  dashed  (E+H  )  curve 
of  Fig.  13  qualitatively  appears  more  reasonable.   For  ex- 
ample, the  peak  in  (E+H  )  at  15°N  appears  to  be  due  to  rela- 
tively strong  evaporative  and  sensible  heat  transfers 
associated  with  the  Trade  Wind  regime,  while  the  higher 
latitude  peak  (35°-50°N)  in  Fig.  13  is  due  to  the  trend  to- 
wards winter  maximum  turbulent-heat  transfers  associated 
with  cool  continental  air  flowing  over  warm  ocean  currents. 
No  attempt  to  select  a  final  K*  value  is  suggested  from 
these  results.   If  the  same  turbulent-transfer  model  can  be 
tested  over  a  large  enough  set  of  seasonal  samples,  it  should 
be  possible  to  assess  the  appropriate  K*-value  which  yields 
the  best  approximation  to  the  atmospheric  heat  balance. 

D.   SURFACE  HEATING  RATES 

Fig.  14  ""depicts  the  zonally-averaged  surface  heating  bud- 
get as  a  function  of  latitude  resulting  from  this  study. 
The  parameters  include  the  net  surface  radiative  warming  (R), 
the  surface  heat  loss  due  to  turbulent  transfer,  -(E+Hr), 
and  the  net  warming  or  cooling  at  the  surface,  (Q ?   +S). 
These  parameters  are  related  by  the  following  expression 

96 


Heating  rafes, 
Surface 


10    15    20   25    30   35   40   45    50   55   60   65 
Latitude  (°N) 
Figure  14.   Surface  heat  budget  disposition.   The  net  warming 

(cooling)  rate  at  the  earth's  surface  (Q   +S)  is  given  by  the 

sum  of  the  net  radiative  surface  warming  rate  R  and  the  loss 

rate  due  to  turbulent  heat  transfer,  -(E+H). 
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LAT 

0 

5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65N 


(a) 

Ri  /  R2 


(b)  (c) 


i  i 

1   2 

vo    1 

VO     2 

.1687/ 

-.0294/ 

.1393/ 

.1720/. 2033 

-.0332/- 

.1539 

.1388/ 

.0494 

.1596/ 

-.0425/ 

.1171/ 

.1463/. 2040 

-.0636/- 

.1786 

.0827/ 

.0304 

.1361/ 

-.0544/ 

.0817/ 

.1276/. 1900 

-.0774/- 

.1748 

.0502/ 

.0152 

.0982/ 

-.1104/ 

-.0122/ 

.0705/. 1482 

-.1701/- 

.1615 

-.0996/- 

.0133 

.0703/ 

-.1544/ 

-.0841/ 

.0499/. 1026 

-.1518/- 

.1254 

-.1019/- 

.0228 

.0335/ 

-.1486/ 

-.1151/ 

.0146/. 0646 

-.0880/- 

.0988 

-.0734/- 

.0342 

.0013/ 

-.0924/ 

-.0911/ 

.0120/ 

-.0188/ 

-.0068/ 

wt. 

Avg 


1037/. 1625 


-.0855/-. 1533 


.0182/  .0092 


Table  V.   Comparisons  of  net  radiative  warming  rate  at  the 
surface  (R),  the  heat  loss  rate  at  the  surface  due  to  turbu- 
lent transport  -(E+H"r),  and  the  net  warming  or  cooling  at  the 

earth's  surface  (Q   +S)  for  16  October  1973  (denoted  by  sub- 

vo 

script  "1").   For  comparison,  annual  mean  values  reported  by 
Malkus  (1962)  are  listed  at  selected  latitudes  (denoted  by 
subscript  "2").   All  values  are  in  ly  min"1. 
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Qvo  +  S  =  R  -  (E+Hr)  (7-4) 

where  S  is  the  oceanic  heat-storage  term.   Again  it  was  not 

possible  here  to  separate  the  contributions,  Q   and  S.   The 

square  points  and  dashed-in  portion  of  the  -(E+IT)  and  Q 

curves  depict  the  effect  of  varying  the  K*  value  in  Eq.  (5-1) 

as  discussed  in  Sec.  VII. C.   The  overall  distribution  in  the 

net  oceanic  heating  rate  (Q   +S)  is  the  same  for  both  choices 

vo 

of  K* .   However  this  rate  is  somewhat  reduced  for  the  case 
of  the  larger  K*  values  suggested  as  applicable  for  the 
tropics . 

Table  V  shows  quantitatively  the  difference  in  the  param- 
eters depicted  in  Fig.  14  as  found  by  this  study  for  16 
October  1973  and  the  climatological  annual  values  reported 
by  Malkus  (1962).   Note  the  smaller  mean  seasonal  radiation- 
al  warming  rate  (R)  and  the  effect  of  the  smaller  mean  heat 
loss  by  turbulent  transfer  in  giving  a  slightly  larger  mean 
warming  of  the  surface. 
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VIII.   CONCLUSIONS 

A  heat  budget  model  for  inclusion  in  the  FNWC  numerical 
weather  prediction  scheme  was  evaluated  using  data  along 
four  meridional  lines  over  the  Pacific  and  Atlantic  Oceans. 
These  data  soundings  were  used  to  compute  a  space-time  aver- 
age of  the  radiational  and  turbulent  heating  parameters  for 
one  day  over  an  ocean  environment.   This  space-time  average 
agreed  reasonably  well  with  values  derived  from  other  stud- 
ies, including  satellite  measurements.   There  were  some  ob- 
vious differences  which  could  be  ascribed  to  the  seasonality 
of  the  data  used  in  this  study  comprising  only  a  one  day 
space-time  average  while  the  comparison  data  was  averaged 
over  years  to  arrive  at  climatological  results. 

Other  differences  could  be  ascribed  to  the  two-cloud 
layer  parameterization  for  cloud  amounts,  the  treatment  of 
cloud  reflectivity  and  absorptivity  for  solar  wave  lengths, 
and  the  parameterization  of  sensible  and  latent  heat  trans- 
fers.  At  this  stage  in  the  model  testing,  it  would  seem 
that  the  radiational  estimates  are  still  somewhat  question- 
able due  to  the  adoption  of  the  Smagorinsky-type  cloud 
amounts.   However  the  regular  variation  of  the  radiative 
terms  from  gridpoint  to  gridpoint  tends  to  support  the  ap- 
parent realism  of  the  radiative  model.   The  turbulent  trans- 
fer heating  rates  of  the  model  seem  to  be  subject  to  a 
greater  degree  of  possible  random  error  due  to  uncertainty 
of  the  vertical  gradients  involved  near  the  surface. 
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To  complete  the  comparisons  of  the  radiational  aspects 
of  this  model  with  cited  climatological  data  (see  Sec.  VII) 
and  with  satellite  data,  further  testing  should  be  done  at 
different  seasons  of  the  year  and  possibly  over  other  oceanic 
areas.   Thus  with  a  full  year  of  representative  data,  refine- 
ments of  the  cloud-cover  parameterization  and  in  the  cloud 
absorption  and  reflection  values  used  in  the  model  may  be 
assessed. 

Finally,  the  greatest  gridpoint  to  gridpoint  variability 
occurred  in  connection  with  the  computation  of  evaporative 
and  sensible  heat  parameters.   It  is  perhaps  here  that  the 
moisture  and  heating  source  terms  used  in  the  FNWC  predic- 
tion model  might  best  be  improved. 
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